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By 
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Chairperson:  C.M.  Van  Vliet 

Major  Department:  Electrical  Engineering 

Noise  properties  of  n"*"  n"  n'*',  p'*'  p"  p'*',  and  p'*’  n"  p'*’  GaAs 
structures  are  investigated.  The  thickness  of  the  n-type  devices  are 
.24um,  .4um,  and  2.0um  while  the  p'*’  n"  p'*'  and  p"*"  p"  p'*'  have  a thickness 
of  .4ym.  These  very  short  devices  are  of  interest  since  they  are  useful 
tools  to  investigate  fundamental  phenomena  such  as  ballistic  transport 
and  1/f  noise. 

Current-voltage  measurements  and  low  frequency  noise  measurements 
at  room  temperature  for  every  device  previously  mentioned  are 
presented.  Comparison  between  the  various  levels  of  1/f  noise  in  the  n- 
type  devices  leads  to  the  conclusion  that  mobility  fluctuations  are 
responsible  for  noise  in  these  devices,  since  in  the  shorter  devices 
(.24um  and  .4um),  the  noise  level  is  low.  This  suggests  some  carriers 
have  a near-ballistic  transport.  The  high  level  of  low  frequency  noise 
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in  the  p'*’  n“  p'*’  device  suggests  the  possibility  of  noise  due  to  traps  in 
this  device. 

The  noise  at  high  frequency  and  at  high  electric  fields  is  further 
investigated  on  the  I.Imiti  n-type  device. 

A step-up  transformer  is  designed  to  bring  the  noise  of  the  device 
to  a level  above  that  of  the  radio  frequency  preamplifier. 

The  noise  at  high  electric  fields  can  be  explained  in  terms  of 
velocity  fluctuations  due  to  intervalley  scattering.  The  electron  drift 
velocity  v(E),  and  the  electron  diffusion  coefficient  D(E)  obtained  from 
the  measurements  by  the  numerical  deconvolution  of  the  I-V 
characteristic  and  the  noise  spectral  density,  are  compared  to  those 
obtained  by  Monte  Carlo  simulation  techniques. 
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CHAPTER  I 
INTRODUCTION 


Progress  in  semiconductor  device  technology  such  as  the  development 
of  the  Molecular  Beam  Epitaxy  (MBE)  made  it  possible  to  fabricate 

devices  with  ever  decreasing  size.  The  use  of  compound  semiconductors 
with  relatively  high  carrier  mobility  such  as  GaAs  to  fabricate  these 
small  devices  has  led  to  the  fabrication  of  low  power  high-speed 
devices.  Beyond  these  practical  considerations  such  small  structures 
are  useful  tools  to  investigate  new  models  of  charge  transport  for  which 
traditional  approaches  no  longer  apply. 

If  one  seeks  to  use  these  small  devices  as  part  of  a detection 

system  like  sensors,  then  the  noise  performance  of  the  device  is  very 
important,  since  the  noise  can  be  a serious  limitation  on  the 
performance  of  the  system. 

Noise  studies  do  not  only  serve  those  practical  purposes.  Low 

frequency  noise  measurements  may  contribute  very  much  to  the 

understanding  of  the  1/f  noise  phenomenon.  Thermal  noise  measurements 
at  higher  frequencies  and  higher  bias  voltages  can  give  useful 
information  about  the  high  electric  field  behavior  of  these  devices. 
All  these  experiments  supported  with  numerical  simulations  are  described 
in  this  work. 

The  organization  of  the  chapters  is  as  follows.  In  Chapter  II  a 
review  of  previous  work  on  noise  measurements  of  submicron  devices  is 
summarized.  Measurement  procedures  of  low  frequency  noise  are  presented 
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in  Chapter  III.  In  Chapter  IV  experimental  results  on  low  frequency 
noise  measurements  on  n'*'  n"  n’’’,  p"'’  p"  p"^,  and  p'*'  n"  p'*’  GaAs  structures 
are  described.  In  Chapter  V an  analysis  and  discussions  of  the  results 
of  the  low  frequency  noise  measurements  are  presented. 

Chapter  VI  is  the  beginning  of  the  second  part  of  this  work,  which 
investigates  the  high  electric  field  noise  properties  of  the  l.lum  n"'"  n“ 
n'*'  GaAs  device.  Measurement  procedures  are  described  in  that  chapter. 
In  Chapter  VII  experimental  results  are  presented.  A Monte  Carlo 
simulation  method  which  is  used  to  calculate  the  electron  drift  velocity 
and  the  electron  diffusion  coefficient  as  a function  of  the  electric 
field  is  described  in  Chapter  VIII.  Chapter  IX  is  devoted  to  the 
discussion  of  the  experimental  results  and  numerical  calculation  of  the 
thermal  noise.  Finally  conclusions  and  recommendations  for  further  work 
are  presented  in  Chapter  X. 


CHAPTER  II 
PREVIOUS  WORK 


In  the  first  part  of  this  chapter  theoretical  work  about  current- 
voltage  characteristics  is  briefly  described.  The  second  part  presents 
the  experimental  work  on  .4ufn  n'*’  n"  n"''  GaAs  and  .47ym  n'*'  p"  n'^  performed 
by  R.R.  Schmidt  [1],  Shur  and  Eastman  [2]  used  the  theory  of  vacuum 
tubes  to  derive  a simple  model.  The  current  density  is  given  by 

J = q n V (2.1) 

where  q is  the  electronic  charge,  n is  the  free  electron  density  and  v 
is  the  velocity. 

The  velocity  v can  be  obtained  from  the  energy  conservation 
equation 


qV^  = 1/2  m*  v^  (2.2) 

where  is  the  potential  of  x due  to  the  applied  voltage  V,  and  m*  the 
effective  mass  of  the  electron. 

The  relation  between  the  potential  and  the  electron  density  is 
given  by  Poisson's  equation.  They  arbitrarily  set  the  initial  velocity 
at  the  cathode  and  the  field  equal  to  zero. 

3/2 

This  approach  leads  to  a 3/2  power  law  (0  ~ V ) for  large  bias 
voltages  (Child's  law)  and  an  1/2  power  law  for  low  bias  voltages. 
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The  same  authors  [3]  later  modified  the  previous  models  by 
introducing  a friction  term  which  takas  into  account  the  few  possible 
scattering  events.  Equation  (2.2)  is  then  replaced  by 

d m*  V 

— (m*  v)  = q E (2.3) 

dt  T 

where  t is  a momentum  relaxation  time. 

Basically  both  models  have  been  questioned  by  the  fact  that  the 
boundary  conditions  were  arbitrarily  set.  Rosenberg  et  al . [4]  have 
shown  that  in  the  case  of  very  short  devices  the  boundary  conditions  are 
very  important  in  modelling  the  devices. 

Holden  and  Debney  [5]  with  their  improved  theory  of  ballistic 
transport  in  one  dimension,  present  an  interesting  approach  to  the 
problem.  They  introduce  the  fact  that  the  electrons  have  an  initial 
velocity  distribution  at  the  contacts.  The  method  of  analysis  of  the 
problem  is  essentially  that  used  for  vacuum  tubes  by  Fry  [6].  This 
method  is  also  known  as  the  Langmuir-Fry  method.  They  recognize  that, 
for  low  bias  at  least,  a potential  minimum  exists  between  the  electrodes 
because  the  cathode  contact  can  supply  more  electrons  to  the  n"  region 
than  reach  the  anode.  An  advantage  of  this  particular  formulation  is 
the  natural  way  in  which  the  boundary  conditions  for  the  problem  are 
imposed.  These  are  specified  by  the  potential  at  the  anode  and  the 
cathode  and  the  total  electron  injection  rates  from  the  n"^  contacts, 
together  with  the  implicit  condition  of  zero  electric  field  at  the 
potential  minimum.  This  situation  is  depicted  in  Figure  2.1.  x‘  is  the 
distance  from  the  cathode  where  the  potential  minimum  occurs  and  V'  is 
its  magnitude. 
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Figure  2.1:  Potential  profile  of 

the  Hoi den-Debney  model 


Anode 


6 


The  problem  is  then  to  describe  the  charge  density  p(x)  in  the 
region  (a)  between  the  cathode  and  the  potential  minimum  and  in  the 
region  (3)  between  the  potential  minimum  and  the  anode. 

For  each  region  Poisson's  equation  is  solved. 

d^V  p(x) 

— 2"  ^ (2.4) 

dx  ee 

0 

In  the  region  (a)  p(x)  is  the  sum  of  four  contributions: 
p^^:  The  density  of  electrons  emitted  by  the  cathode  which  have  enough 

energy  to  overcome  the  potential  minimium  and  reach  the  anode.  They 
have  a velocity  greater  than  a critical  velocity  vq  such  that 


1/2 

”*  'o 

- qV  > 

0 

(2.5) 

or 

< 

O 

V 

2q 

(— V 
m* 

■)l/2  ^ 

''o 

(2.6) 

CO 

n (v  ) 

then 

= 

q / 

0 

CL 

< 

o 

(2.7) 

''o 

V 

The  density  of  electrons 

emitted 

by  the 

cathode  which  do  not  have 

enough  energy  to  pass  the  potential  minimum.  Thus,  there  are  two 
equally  dense  streams  of  electrons  passing  in  opposite  directions  and 
having  kinetic  energies  less  tha  qV  or  in  other  words 


0 < Vq  < Vg  (2.8) 


then 
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Po  = 2q  fr 


— 0 


(2.9) 


P3:  The  density  of  electrons  emitted  from  the  anode  which  have  enough 

energy  to  overcome  the  potential  minimum  v‘ 


pQ  = q / 


P(Vn) 


dv. 


'^0  '' 


(2.10) 


where 


2q 

v"  = [_  (V,  - V')] 


1/2 


0 


m* 


Vj,:  Anode  voltage. 

p^:  The  background  space  charge  due  to  ionized  donors 


P4  = 


N 


D 


(2.11) 


(2.12) 


For  the  region  (8)  the  same  expressions  hold  for  p except  for  p2  which 
becomes 


n(vn) 


= 2q  / 


dv 


0 

1/2 


(2.13) 


P2:  The  density  of  the  electrons  emitted  from  the  anode  which  have 

enough  energy  to  escape  from  the  anode  but  not  enough  to  pass  the 
potential  minimum. 
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In  the  previous  expressions  of  p,  q is  taken  as  equal  to 
-19 

-1.6  X 10  C . In  each  case  Poisson's  equation  is  solved  by  multi- 
plying both  sides  of  equation  (2.4)  by  2 ^ and  integrating. 

We  have 


dV  p 2 

(— ) = / P^-  dV  (2.14) 

dx  eGq 

where  has  been  described  before.  In  addition  to  Poisson's  equation, 
the  current  density  equation  is  needed. 

.r  !a 

j = q N e ''t  (1  - e '^T)  (2.15) 

+ UT 

where  N = Ng  Vg  and  Vj  is  the  thermal  voltage,  — . These  equations 
must  be  solved  sel f-consi stently  to  obtain  the  anode  voltage  for  any 
specified  net  current  density  j flowing  between  the  cathode  and  the 
anode.  According  to  Holden  and  Debney,  it  is  difficult  to  identify  the 
ballistic  motion  of  the  electrons  by  the  current-vol tage  character! Stic 
only.  R.R.  Schmidt  [1]  measured  the  current  as  a function  of  the  bias 
voltage  and  the  low  frequency  noise  for  different  bias  currents  on  .4um 
n"*"  n“  n"^  and  .47ym  n'*'  p“  n'*’  devices. 

The  I-V  characteristic  of  the  n-type  is  linear  up  to  about  1 volt 
for  300  and  77K.  Neither  the  nor  the  dependence  was 

observed.  He  concluded  that  the  theory  of  Shur  and  Eastman  [3]  is  not 
adequate  to  describe  his  results.  He  believed  that  a more  realistic 
model  is  that  of  Holden  and  Debney  which  gives  a dependence  at 

high-bias  current  for  a .Sum  device.  The  low  bias  voltage  regime  can  be 
explained  by  the  thermionic  and  diffusion  drift  model  of  van  der  Ziel  et 
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a1 . [7].  In  that  model  they  calculated  separately  the  resistance  due  to 
the  thermionic  emission  and  due  to  the  diffusion  and  drift. 

For  the  1/f  noise  R.R.  Schmidt  [1]  reported  a Hooge  parameter  a as 

O 

low  as  5.10"°  at  room  temperature  for  the  .4um  n-type  device.  This 
value  is  five  orders  of  magnitude  smaller  than  the  value  of  6 x 10"^ 
reported  by  Hooge  et  al . [8]  for  n-type  bulk  GaAs.  The  noise 

measurement  for  different  temperatures  (77K  and  300K)  suggested  that  the 
device  is  not  purely  ballistic,  but  the  very  low  value  of  ct  suggests 
near-ballistic  transport.  The  current-voltage  characteristic  of  the  n"^ 
p"  n'*’  device  is  much  more  complicated  than  that  of  the  n'*’  n"  n'*' 
device.  It  shows  substantial  nonlinearity.  There  are  two  regimes  with 
a transition  region  in  between.  At  low  bias,  the  characteristic  is 
linear  with  a DC  conductance  which  is  temperature  dependent  and  100 
times  smaller  than  that  of  the  n-type  device.  The  high  bias  regime  is 
temperature  independent  and  is  also  linear  with  a larger  conductance 
than  the  low  bias  regime.  R.R.  Schmidt  [1]  explained  the  low  bias 
regime  as  follows.  The  conductance  depends  on  temperature  because  the 
hole  concentration  does.  At  high  bias  the  transport  is  near-ballistic 
and  thus  temperature  independent  as  observed  in  the  n-type  device.  He 
explained  the  difference  of  conductance  in  the  high  bias  regime  (the  p- 
type  device  has  a conductance  four  times  smaller  than  that  of  the  n-type 
device)  by  the  fact  that  in  the  p-type  case,  the  barrier  is  higher.  The 

noise  picture  fits  well  in  the  above  description.  At  low  bias  R.R. 

2 

Schmidt  [1]  measured  a noisiness  factor  S^j(f)f/I  which  was  4 x 10° 
times  that  for  the  n-type  device  value.  At  high  bias  the  noise  is  much 
smaller.  This  was  shown  by  measurements  performed  by  C.F.  Whiteside 
[9].  The  low  value  of  the  noise  corresponds  to  the  near  ballistic 


10 


transport  suggested  for  high  bias.  In  Chapter  III,  the  measurement 
procedures  of  1/f  noise  for  low  impedance  devices  are  described. 


CHAPTER  III 

MEASUREMENT  PROCEDURES 


3.1  Current  Voltage  Measurements 


In  order  to  measure  the  DC  I-V  characteri stic  of  each  device, 
pulsed  measurements  were  done  to  prevent  excess  heating  which  may  change 
and  ultimately  damage  the  devices.  An  HP  214A  pulse  generator  was  used 
to  produce  pulses  of  lus  width  with  a repetition  rate  of  lOOHz.  The 
set-up  is  shown  in  Figure  3.1.  The  current  flowing  through  the  device 


I 


(3.1) 


where  R^  is  a biasing  resistance  whose  value  should  be  chosen  carefully 
so  that  reasonable  voltages  can  be  measured  across  the  device. 


3.2  Noise  Measurements 


3. 2. a Method 

The  noise  measurement  procedure  is  the  usual  three  points 
measurement: 

A D.C.  voltage  is  applied  to  the  terminals  of  the  device 
under  test. 

A noise  calibration  source  is  applied  to  the  terminals 
of  the  small  signal  equivalent  resistance  of  the  device. 

No  bias  voltage  is  applied  to  the  small  signal 
equivalent  resistance  of  the  device. 
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Figure  3.1:  Current-vol tage 

measurement  set-up 
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The  fourth  information  is  the  magnitude  of  the  calibration  itself.  This 

information  is  sufficient  to  compute  the  excess  noise  of  the  device. 

This  is  easily  understood  by  analyzing  the  following  block 

diagram.  In  position  1 the  device  is  connected  to  a battery  through  a 

calibration  resistance  R3.  In  term  of  noise,  the  battery  does  not 

contribute  at  all.  Its  Norton  equivalent  is  an  open  circuit.  The 

device  under  test  can  be  represented  by  a current  source  expressing  the 

excess  noise  in  parallel  with  the  small  signal  resistance  of  the  device 

itself.  The  low  noise  amplifier  is  modelled  as  a noiseless  amplifier 

with  a voltage  gain,  a voltage  noise  source  in  series  at  the  input  and 

the  input  impedance  and  a current  noise  source  in  parallel.  Finally, 

the  calibration  resistance  R^,  the  input  impedance  of  the  amplifier  and 

the  device  resistance  can  be  lumped  together  into  an  equivalent 

resistance  R^^  = R^^^  ||  R^^  ||  R^  representing  an  equivalent  current 
4k  T 

noise  source  . The  meter  reading  can  be  expressed  in  terms  of  the 

eq 

power  spectral  densities  of  the  noise  sources  as  follows; 


m2  = 


B [S 


Vamp 


+ (S 


4k  T 


Idev 


Sy  ) 

lamp' 


R^  ] 
eq^ 


(3.2) 


eq 


where  B is  the  bandwidth  of  the  spectrum  analyzer.  In  position  2 the 
device  is  connected  to  a calibrated  voltage  noise  source  through  the 
same  source  resistance  R5.  The  Norton  equivalent  is  a current  source  in 
parallel  with  R5.  The  device  under  test  is  replaced  by  a resistor  whose 
resistance  is  equal  to  the  small  signal  resistance  of  the  device.  This 
is  very  important  for  the  case  of  nonlinear  devices. 

The  expression  of  the  meter  reading  follows: 
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B [S 


Vamp 


+ (- 


veal 


4k  T 

+ + S 

’^eq 


lamp'  ^eq-* 


(3.3) 


Finally  in  position  3 the  calibration  resistance  Rj  is  grounded  and  the 
device  is  replaced  by  its  small  signal  equivalent.  The  meter  reading  is 
then  given  by 


4k  T 


+ S,  ) R^  ] 
lamp  eq 


(3.4) 


eq 


Equations  3.2  and  3.4  yield 


B R^  S.  . 

13  eq  Idev 


(3.5) 


similarly  3.4  and  3.4 


M 


2 

2 


^Vcal 


(3.6) 


From  equations  3.5  and  3.6  one  obtains  the  excess  noise  of  the  device 


Idev 


^Vcal 


(3.7) 


Every  quantity  of  the  right  hand  side  of  equation  3.7  is  known.  This 
method  determines  in  a unique  way  the  excess  noise  of  the  device  under 
test. 


3.2.b  Description  of  the  Set-Up 

In  this  paragraph,  the  equipment  shown  in  the  block  diagram  of 
Figure  3.2  will  be  described.  A Hewlett  Packard  Model  3582A  is  used  as 
wave  analyzer.  This  model  is  a dual-channel  spectrum  analyzer  covering 
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Low  Noise 
Amp  I if  ier 


Figure  3.2:  Block  diagram  for  the 

noise  measurement  set-up 
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the  frequency  range  .02Hz  to  25KHz,  featuring  fast  Fourier  transform  for 
both  channels.  Its  capability  of  measuring  input  signals  from  +30db 
down  to  -120db  is  very  important  when  one  deals  with  small  signals  like 
noi se. 

The  low  noise  preamplifier  used  is  the  one  called  PAl  designed  by 
R.R.  Schmidt  [1].  The  circuit  diagram  of  this  amplifier  is  shown  in 
Figure  3.3.  Details  of  this  configuration  are  well  described  in  R.R. 
Schmidt's  Ph.D.  thesis  [1].  The  main  features  that  make  this  amplifier 
very  useful  for  low  frequency  and  low  impedance  device  noise  measurement 
are  its  extended  low  frequency  response,  its  high  voltage  gain  of  the 
order  of  90db  and  its  relative  low  input  impedance.  This  low  input 
impedance  results  from  the  fact  that  p n p transistors  are  used  which 
have  low  base  resistances.  The  calibration  source  is  provided  by  the 
spectrum  analyzer  itself.  It  produces  a computer  generated  pseudo 
random  signal.  This  pseudo  random  signal  has  a flat  frequency  response 
over  the  range  of  the  frequency  of  interest. 

In  Chapter  IV  results  on  current  voltage  measurements  and  low 
frequency  noise  measurements  are  presented  for  n'*’  n"  n'*',  p"*"  n"  p"^,  and 
p'*'  p“  p'*’  GaAs  structures. 
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Figure  3.3: 


Low  noise  preamplifier 


CHAPTER  IV 
MEASUREMENT  RESULTS 


4.1  Current-Voltage  Characteristic 
4.1. a The  n'*'  n~  n'*'  Device 

The  diodes  are  mesa  structures  with  a lightly  doped  n-type  layer 
sandwiched  between  heavily  doped  n'*’  layers  which  serve  as  contacts.  The 
structure  is  shown  in  Figure  4.1.  The  n'*’  regions  have  a doping  density 
of  lO^''^  cm"^  while  the  active  region  has  a doping  density  of  2 x 10^^ 
cm”'^.  The  diameter  of  the  mesa  is  approximately  lOOym.  These  n-type 
devices,  as  well  as  the  others  used  in  the  measurements,  were 
manufactured  by  molecular  beam  epitaxy  at  the  Cornell  University 
Submicron  Research  and  Resource  Facility.  The  measurements  were 
performed  on  devices  with  .24um,  l.lym,  and  2.0ym  thickness.  The  DC 
current-voltage  characteristic  of  the  l.lum  n'*’  n"  n'*'  device  at  room 
temperature  is  shown  in  Figure  4.2. 

At  low  bias,  the  characteristic  is  linear.  The  device  behaves  like 
a resistor.  At  high  bias,  around  .2V,  the  conductivity  starts  to 
decrease.  This  corresponds  to  the  transfer  of  electrons  to  the  upper 
valley  of  the  conduction  band  where  the  electrons  have  a large  effective 
mass  and  thus  a lower  mobility.  The  l.lym  device  behaves  like  a bulk  n- 
type  GaAs  semiconductor.  The  2.0iiin  device  has  a slightly  different 
doping  profile.  The  n’'*  regions  have  a doping  density  of  2.10^®  cm“^. 
The  DC  current-vol tage  characteristic  at  room  temperature  is  shown  in 
Figure  4.3.  The  characteristic  is  linear  at  bias  voltages  below  .2V. 
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100  JLim  diameter 


Figure  4.1:  l.lym  n-type  device 
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Figure  4.2:  Current-voltage  characteristic  of  the 

l.lym  n-type  device  at  room  temperature 
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Figure  4.3:  The  current-voltage  characteristic  of 

the  2.0  ym  n-type  device  at  room  temperature 
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Figure  4.4:  Current-voltage  characteristic  of  the 

.24ym  n-type  device  at  room  temperature 
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For  higher  bias  voltages  the  characteristic  is  quadratic,  the  decrease 
of  conductivity  as  observed  in  the  l.luin  device  was  overshadowed  by  the 
effect  of  space  charge  injection. 

The  I.V  characteristic  of  the  .24um  device  is  shown  in  Figure 
4.4.  The  regime  is  linear  for  both  low  and  high  bias  voltage.  This  is 
similar  to  the  .4um  device  measured  by  Schmidt  [1].  This  may  indicate  a 
near  ballistic  transport. 

4.1. b  The  p'*'  p~  p'*'  Device 

The  p-type  device  has  the  following  structure:  a .55um  p'’”''  region 
of  10^^  cm"^  doping  density  is  grown  on  a p'*’  substrate  followed  by  a 
.55um  p 10^°  doped  region,  a .4ym  lightly  doped  p-region  and  finally  a 
.55um  p'*’  10^^  and  a .55ym  p''”''  10^^  to  serve  as  contact  on  top  of  the 
mesa.  It  is  really  a p"'"'’  p'*'  p"  p'’’  p''"*'  structure  rather  than  a p'*'  p“ 
p'*’.  The  DC  current-vol tage  characteristic  at  room  temperature  is  shown 
in  Figure  4.5.  The  characteristic  is  linear  with  a resistance  of 
1.5^2.  The  device  behaves  like  a pure  resistance. 

4.1. C  The  p~*~  n~  p"*"  Device 

The  structure  of  the  p'*'  n"  p'*’  device  is  similar  to  that  of  the  p'^ 
p“  p"^  device  except  that  the  region  .4ym  p-region  is  now  replaced  by  a 
.4ym  n-region  with  a doping  density  of  2 x 10^^  cm~^.  The  dopent  is 
silicon.  The  DC  current-vol tage  character! Stic  at  room  temperature  is 
shown  in  Figure  4.6.  The  characteristic  is  linear  up  for  bias  voltages 
up  to  .05V  then  a 3/2  power  low  regime  is  observed  up  to  IV.  Beyond  IV 
it  seems  that  the  characteristic  is  linear  again.  However,  it  is 
difficult  to  really  assume  that  just  with  two  experimental  data 
points.  It  was  impossible  to  go  beyond  2V  with  the  available  equipment 
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Figure  4.5:  The  current-vol tage  characteristic  of 

the  .4ym  p-type  device  at  room  temperature 


25 


Figure  4.6:  The  current-voltage  characteristic  of 

the  ,4  P+  n"  p+  device  at  room  temperature 
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since  the  pulse  generator  requires  a high  impedance  load  to  be  able  to 
produce  undistorted  output  signals. 

4.2  Low  Frequency  Noise 

4. 2.  a n'*'  - n~  n'*'  Device 

The  low  frequency  excess  noise  for  the  device  was  measured 

from  1 to  25KHz.  The  current  noise  spectral  densities  for  several 
bias  currents  at  300K  is  shown  in  Figure  4.7.  The  frequency  dependence 
of  the  excess  noise  is  1/f.  From  the  different  bias  current  levels,  the 
current  noise  spectral  density  Sj  is  proportional  to  the  square  of  the 
bias  current.  This  is  the  expected  behavior  of  the  1/f  noise  for  a 
device  with  a linear  current-voltage  characteristic.  The  current  noise 
spectral  density  at  room  temperature  of  the  2.0ym  device  for  several 
bias  currents  is  plotted  in  Figure  4.8.  The  frequency  dependence  is 
again  1/f  at  low  frequencies  up  to  IKHz.  A steeper  slope  is  then 
observed  beyond  IKHz.  The  current  noise  spectral  density  is 

proportional  to  the  square  of  the  bias  current  as  expected  since  the 
magnitude  of  the  bias  currents  are  still  within  the  linear  part  of  the 
current-voltage  characteristic. 

The  current  noise  spectral  density  at  room  temperature  for  the 

.24um  device  is  shown  in  Figure  4.9. 

4.2. b  The  p~*~  p~  p'^  Device 

The  low  frequency  excess  noise  of  the  .4ym  p-type  device  was 

measured  from  IHz  up  to  25KHz  at  room  temperature.  Figure  4.10  shows 
the  current  noise  spectral  density  as  a function  of  frequency.  The 
magnitudes  of  the  noise  spectral  density  decreases  inversely  with  the 

frequency.  The  spectral  density  is  again  proportional  to  the  square  of 
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Figure  4.7:  Low  frequency  current  noise  spectra  for  the 

l.lym  n-type  device  at  room  temperature 
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Figure  4.8:  Current  noise  spectra  of  the  Z.Oym 

n-type  device  at  room  temperature 
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Figure  4.9:  Current  noise  spectra  of  the  .24ym 

n-type  device  at  room  temperature 
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Figure  4.10;  Current  noise  spectra  of  the  .4ym 

p+  p-  p+  device  at  room  temperature 
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the  bias  current  as  expected  since  the  current-vol tage  characteristic  of 
the  device  is  linear. 

4.2.C  The  p'*'  n~  p~*~  Device 

The  low  frequency  noise  of  the  .4um  p'*’  n“  p'*'  device  was  measured  in 
the  same  range  of  frequencies  as  previously.  Figure  4.11  shows  the 
noise  spectral  density  as  a function  of  the  frequency  at  room 
temperature  for  different  bias  currents.  The  spectra  go  as  1/f  at  least 
for  low  biases.  For  low  bias  currents,  the  noise  is  proportional  to  the 
square  of  the  bias  current,  but  as  the  bias  current  increases,  this  is 
not  true  anymore  since  the  current-voltage  character! stic  is  not  linear. 

In  the  next  chapter  these  experimental  results  will  be  analyzed  and 
compared  with  existing  theories. 


ca2/hz:> 
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Figure  4.11:  Current  noise  spectra  of  the  .4ym 

p+  n-  p+  device  at  room  temperature 


CHAPTER  V 

HOOGE'S  PARAMETER  AND  DISCUSSION  OF  THE  RESULTS 


5.1  Quantuin  1/f  Noise  Theory 

As  a starting  point  we  use  Handel's  formula  [10]  which  expresses 
the  relative  spectral  density  of  the  scattering  cross  section  as 


Sa 


<a> 


2 


2cxA 


f 


(5.1) 


for  low  frequency  photon  emission,  one  obtains 

2a  Av  y 

ccA=  — (— )^  (5.2) 

3ir  C 


where  a=  1/137  is  the  fine  structure  coefficient.  If  the  current 
fluctuation  is  solely  due  to  mobility  fluctuation  then  one  can  write 


Sj(f) 


(5.3) 


According  to  Hooge  [11] 

Si(f)  _ 

Nf 


(5.4) 


where  f is  the  frequency,  I the  DC  bias  current,  N the  number  of 
carriers,  and  a,^  the  Hooge's  parameter. 
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Several  different  scattering  mechanisms  determine  the  net  mobility 
of  an  electron  and  it  is  given  by 


1 1 
y i 


(5.5) 


or  in  terms  of  spectral  density 


(5.6) 


From  equations  (5.3)  and  (5.4)  and  for  one  type  of  scattering 
mechani sm 


The  overall  aj_j  is  then 


(5.7) 


^ "‘i  (5.8) 

i y . 

1 

Kousik  [12]  has  derived  expressions  of  5^  /y^. ^ for  different 

scattering  mechanisms,  such  as  impurity  scattering  polar  optical  phonon 
scattering,  and  acoustical  phonon  scattering  and  consequently  the 
individual  a.,  by  relating  the  mobility  fluctuations  to  the 

fluctuation  of  the  scattering  cross  section  given  by  Handel's  formula. 
This  was  achieved  by  considering  the  fluctuations  of  the  relaxation  time 
associated  with  each  collision  process. 
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5.2  Experimental  Determination  of  for  N-Type  Device 

The  next  step  is  now  to  extract  this  Hooge's  parameter  from  the 
experimental  measurements.  Using  again  Hooge's  empirical  formula,  we 
have 


a 


H 


Si(f) 


(5.9) 


We  need  to  calculate  correctly  N.  Van  der  Ziel  and  Van  Vliet  [13]  gave 
a formula  to  calculate  N,  which  is  also  valid  for  non-homogeneous 
samples  or  mesa  structures 


Sj(f) 


1 1 I dx 

Ji 

f AL  L ^ n(x) 


(5.10) 


where  L is  the  total  length  of  the  device,  A the  cross  sectional  area, 
and  n(x)  the  electron  density  in  the  active  region.  The  problem  has 
always  been  the  calculation  of  the  integral  since  one  does  not  know 
exactly  n(x).  Schmidt  [1]  used  an  average  value  n which  he  assumed 
equal  to  twice  the  value  of  the  donor  density  Nq. 

For  a device  with  m layers  we  have 


L = m L' 


(5.11) 


where  L'  is  the  length  of  one  layer. 

Substituting  in  equation  (5.11)  we  have 


Sj(f) 


1 1 


rL' 


dx 


fm  AL‘  L 


n(x) 


(5.12) 
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or 


Sj(f) 


f m N 


eff 


(5.13) 


where  we  defined  as 


1 11  I I dx 

= /q  

AL*  L‘  ^ n(x) 


(5.14) 


Instead  of  approximating  n(x)  by  an  average  value,  we  computed 
numerically  the  carrier  concentration  profile  n(x)  and  evaluated  the 
integral  in  equation  (5.14).  That  was  done  by  simultaneously  solving 
Poisson's  equation  for  the  entire  device  and  the  current  continuity 
equation.  The  program  was  developed  by  S.  Tehran!  [14].  These  carrier 
concentration  profiles  are  shown  in  Figure  5.1,  5.2,  and  5.3  for  the 
.24,  .4,  and  1.1  ym  devices  respectively.  The  effect  of  carrier 

spillover  is  important  in  the  short  devices.  For  the  .24  ym  device,  the 
carrier  concentration  minimum  is  around  10^®  cm"^  which  is  ten  times  the 
background  doping  density.  The  .4  ym  device  showed  a minimum  of  3.10^^ 
cm"'*  which  is  three  times  the  background  doping  density.  Finally  the 
1.1  ym  device  is  not  affected  by  carrier  spillover. 

In  Table  5.1  computed  according  to  equation  (5.13)  is  listed 

for  the  three  devices. 


Table  5.1:  Hooge's  parameters  for  the  n-type  devices 


Length 
.24  ym 
.4  ym 


^eff 


1. 

10-15 

3.4 

10^ 

1.7 

10-^ 

1.7 

10-15 

1.73 

10^ 

1.5 

10-7 

2.6 

10-14 

2.4 

10^ 

6.2 

10-7 

m 

5 

5 


1.1  ym 


1 
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Figure  5.1:  Carrier  concentration  profile 

for  a .24um  n-type  device 
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Figure  5.2:  The  carrier  concentration  profile 

for  the  .4um  n-type  device 
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Figure  5.3:  Carrier  concentration  profile  of 

the  l.lym  n-type  device 
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The  low  values  of  aj^  for  the  .24vim  and  the  .4um  devices  indicate 
the  reduction  of  noise  due  to  near  ballistic  transport. 

Kousik  [12]  performed  some  preliminary  calculations  on  each 

-8 

individual  a.  and  using  equation  (5.8)  he  found  = 2.10"  for  n 
type  GaAs  at  room  temperature,  considering  the  impurity  scattering, 
acoustical  phonon  scattering,  and  polar  optical  scattering.  Comparing 

this  value  to  those  from  Table  1,  they  differ  by  a factor  of  ten.  This 
can  be  explained  by  some  uncertainty  in  the  calculation  of  Ng^f,  the 

uncertainty  inherent  to  any  measurement  or  the  uncertainty  in  the 
material  parameters  used  to  compute  the  theoretical  a. 

These  experimental  values  and  the  theoretical  one  computed  by 
Kousik  [12]  make  us  strongly  believe  that  in  this  case  the  mobility 

fluctuation  is  responsible  for  the  1/f  noise  in  these  samples  and  that 
can  be  understood  by  using  Handel's  theory  [10].  The  same  calculation 
has  been  carried  out  for  the  2.0  um  n-type  device. 

The  value  obtained  for  is 

= 6.2  10"^  (5.15) 

This  value  is  three  orders  of  magnitude  larger  than  the  theoretical 
value.  This  can  not  be  explained  by  the  quantum  1/f  noise  theory.  This 
device  may  have  a different  source  of  noise  other  than  the  fundamental 
noise  due  to  mobility  fluctuation.  A closer  look  to  the  spectrum  may 
suggest  the  possibility  of  generation-recombination  noise  since  the 
spectrum  show  1/f^  slope  beyond  1 kHz. 
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5.3  The  p~  p'*'  Device 

From  the  data,  we  obtained  a value  of  an 

H 

a^  = 2.6  10'^  (5.16) 

for  the  .4  ym  p-type  device  at  room  temperature.  Even  though  the 

theoretical  value  of  for  a p-type  GaAs  material  has  not  yet  been 

calculated,  this  experimental  value  is  not  likely  to  be  explained  by  the 
quantum  1/f  noise  theory.  One  has  to  consider  other  possible  sources  of 
noise.  Contacts  on  p-type  epi layer  are  not  as  good  as  those  of  teh  n- 
type  one,  and  these  may  give  rise  to  a source  of  1/f  noise. 

5.4  The  p~^  n~  p'*'  Device 
The  value  of  aj_|  in  this  case  is 

a^  = 2.5  10'^  (5.17) 

which  is  again  too  high  to  be  explained  by  the  quantum  1/f  noise 
theory.  This  value  is  even  ten  times  larger  than  the  value  for  the  p'*’ 
p"  p'*'  device  which  is  thought  to  produce  contact  noise.  A closer 

comparison  with  the  previous  work  done  by  Schmidt  [1],  however,  reveals 
that  the  noise  has  a similar  magnitude  to  that  of  the  0.47  ym  n'*'  p"  n'*’ 
device.  Schmidt  [1]  did  not  calculate  a.^  , rather  he  used  the  quantity 
(Sj/I  )f  . An  assumption  that  the  noise  behaves  like  1/f  would  give  an 
of  1.4  10“^. 

The  n'*’  p“  n'*’  and  the  p"^  n~  p'*’  devices  are  not  only  similar  with 
respect  to  low  frequency  noise,  but  also  with  respect  to  current-voltage 
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characteristic.  Both  exhibit  a linear  regime  at  low  voltages,  a faster 
increase  in  conductivity  and  another  linear  regime  beyond  IV.  That 
second  linear  regime  is  not  obvious  though  for  the  .4  ym  p'*’  n" 
device.  As  in  the  n'*’  p"  n'*’  case,  we  can  explain  the  noise  in  terms  of 
trap  filling  noise.  This  type  of  noise  would  easily  overshadow  true  1/f 
noi se. 


5.5  Conclusion 

We  believe  that  the  results  from  Section  5.2  are  a real  step  toward 
the  understanding  of  the  1/f  noise  phenomenon.  In  that  section  we 
established  two  major  facts:  first,  the  Hooge's  parameter  does  not 

really  depend  on  the  length  of  the  device  for  devices  with  L < 1 m. 
Second,  for  these  devices,  the  experimental  values  for  are 

comparable  to  the  theoretical  ones  based  on  the  theory  of  quantum  1/f 
noise.  There  is  one  order  of  magnitude  difference  but  that  is  easily 
accounted  for  by  the  uncertainties  we  have  already  discussed.  In  that 
respect,  the  measurements  have  been  conclusive  to  the  understanding  of 
the  fundamental  1/f  phenomena.  This  work  really  establishes  the  idea 

that  1/f  noise  in  semiconductors  can  be  explained  by  mobility 
fluctuation  using  Handel's  theory  [10]  as  elaborated  by  Kousik  and  Van 
Vliet  [12].  The  lack  of  success  with  the  p'*’  p"  p'^  device  was  explained 
in  Section  5.3.  The  2.0  ym  n-type  device,  the  p“  n"^,  and  the  p"^  n" 
p"'’  devices  have  different  noise  mechanisms  in  addition  to  1/f  noise. 


CHAPTER  VI 

MEASUREMENT  PROCEDURES  FOR  PULSED  NOISE  MEASUREMENTS 


The  noise  measurement  in  this  part  is  completely  different  from  the 
method  used  in  the  first  part  of  this  work.  A general  description  of 
the  method  is  given  in  the  first  paragraph  then  a detailed  description 
of  the  equipment  used  is  given  thereafter. 

6.1  Experimental  Method 

The  block  diagram  of  the  experimental  set-up  is  shown  in  Figure 

6.1.  A pulsed  bias  method  was  used  to  avoid  excessive  Joule  heating  and 
possible  destruction  of  the  device  at  high  electric  fields.  The  pulse 
generator  is  synchronized  with  a switch  at  the  input  of  the  power 
detector  so  that  one  only  measures  the  noise  when  the  bias  is  applied. 
Bareikis  et  al . [15]  have  used  this  method  before,  however,  a detailed 
description  of  the  step-up  transformer  method,  suggested  by  Dr.  van  der 
Ziel,  at  the  input  of  the  amplifier  is  given  here.  This  network  is  very 
important  because  the  device  has  a very  low  impedance  (<ln)  compared  to 
the  50Q  input  impedance  of  the  amplifier.  As  a consequence,  the  device 
noise  would  drown  in  the  amplifier  noise  without  a step-up  network. 

6.1.  a Step-Up  Transformer 

Since  the  amplifier  has  a noise  equivalent  resistance  of  60n  and 
the  device,  a resistance  less  than  Ifi,  a tuned  circuit  is  needed  to 
increase  the  resistance  seen  by  the  amplifier.  The  noise  equivalent 
circuit  is  shown  in  Figure  6.2,  where  r is  the  series  resistance  of  the 
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Generator  
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Figure  6.1:  Block  diagram  of  the  experimental  set-up 
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Figure  6.2:  Noise  equivalent  circuit  of  the 

step-up  transformer 
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inductor,  r^  the  device  small  signal  resistance,  and  v/4kTrAf  and 
/4kTr^Af  (k  is  the  Boltzmann  constant)  are  the  thermal  noise  sources 
associated  with  r and  r^  respectively. 

The  impedance  at  the  terminals  is  given  by 

juL  + r + r^ 

2 (6-1) 

1 - CO  CL  + jcuC  (i"  + r^) 

If  one  operates  at  the  tuning  frequency  cog  = and  if  one  chooses  (oq 
so  that  r + r^  <<  |ju)gL| , then 


L 

z « _ 

C 


1 

(r  + r^) 


(6.2) 


The  spectral  density  S of  the  AC  open-circuited  voltage  fluctuations  at 
the  output  of  the  network  is 


S = 4k(Tgr  + T^rJ 


jUgC 


d d' 


“i — r + r + r . + jo)  L 

jcJnC  d '^0 


(6.3) 


where  Tg  is  the  inductor  temperature  and  the  device  temperature. 
Carrying  out  the  algebra  and  again  using  the  fact  that  cogCL  = 1, 
equation  (6.3)  becomes 


S = 4k  (T„r  ^ 1 r ) ^ (6.4) 

c (r  + r^) 

The  condition  for  "up"  transformation  of  the  noise  equivalent  resistance 
of  the  device  is  easily  derived  from  equation  (6.4)  if  one  neglects  r. 
Then  S becomes 
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S = 4KT. 


Cr 


(6.5) 


'Up"  transformation  takes  place  if 


L L 

> r , or  — > rj 

Cr,  " C “ 


(6.6) 


For  the  series  resonance  step-up  transformer,  it  holds  that 


2 2 
C ^ 


(6.7) 


where  Q is  the  quality  factor  of  the  transformer.  Hence,  condition 
(6.6)  is  fulfilled  as  long  as  Q > 1.  However,  a high  Q circuit  has  a 
narrow  passband;  therefore  transformers  for  each  frequency  of  interest 
(11  and  17MHz)  were  built. 


6.1.b  Determination  of  the  Series  Resistance  of  the  Inductor 
If  one  replaces  the  device  by  a short  circuit,  one  obtains 


S 


1 


(6.8) 


Next,  a known  resistance  (rg)  is  used.  Equation  (6.4)  gives 


S 


2 


L 

C(r  + Cq) 


Equations  (6.8)  and  (6.9)  yield 


(6.9) 


(6.10) 


From  equation  (6.10)  one  obtains  the  value  of  the  series  resistance  r of 
the  inductor  at  the  tuning  frequency. 
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6,l.c  Determination  of  the  Device  Resistance  with  Noise  Measurement 

The  actual  device  is  then  placed  in  series  with  the  inductor.  In 
thermal  equilibrium  with  both  the  device  and  the  inductor  at  room 
temperature  Tq,  one  finds 


L 1 

^3  = k T 

° C (r  + r^) 

Equations  (6.9)  and  (6.11)  give 

^3  ^ + '‘O 

From  equation  (6.12)  one  can  calculate  the  value  of  r^j. 


(6.11) 


(6.12) 


6.1.d  Noise  Measurement  with  Bias  Applied 

All  previous  calculations  are  valid  under  thermal  equilibrium 

conditions.  However,  one  is  interested  in  the  high  field  behavior  of 

the  noise  so  one  needs  to  include  a biasing  circuit  in  the  previous 

model.  The  circuit  is  shown  in  Figure  6.3  and  the  noise  equivalent 

circuit  circuit  in  Figure  6.4  Let  = /4KT70f  , V.  = /4KT  ,r  ,Af  , 

s usd  d d * 

and  = i/4KTQrAf  . Possible  excess  noise  of  the  device  is  taken  into 
account  by  using  the  noise  temperature  T^  in  the  expression  for  V^j.  A 
simple  analysis  of  the  circuit  gives  for  the  voltage  fluctuations  at  the 
AC  open-circuited  output  terminals  at  resonance  frequency  ug 


>„c(Rsi-  + t rr^) 


V 


(6.13) 
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Figure  6.3:  The  step-up  transformer  with  the 

OUT  including  the  pulse  bias  network 
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Figure  6.4:  Noise  equivalent  circuit  of  the  step-up 

transformer  with  the  OUT  including  the 
bias  network 


51 


Since  the  noise  sources  are  not  correlated,  the  mean  cross-products 
^di^s'  '^d'^r  ''^'^''sh.  Hence  one  obtains 


L ''d^^V  '*■  ^ 

S(.J=-  " ^ 


L «s<Sr,  " 5^  ) 
+ _ dr 


C (R^r  t R^rj  + C (R^r  + R^r_j  t r^r)^ 


(6.14) 


L 

+ — 


2 f'd'^s^V 


C (R^r  + R^r^  + r^r) 


Since  r^j  <<  R^  and  r <<  R^  (R^  = 50n),  terms  with  r^  and  r^R^  in  the 
numerator  are  dropped.  Equation  (6.14)  becomes 


L 1 Sy 

5(“o' r"'  r2 

Cr2(lr^)2 
d 


(6.15) 


with  Sy  = 4kTgr  and  Sy  = 4kT^r^  . Equation  (6.15)  can  be  written  as 


S(o)o)  = 4k 


r,Z  ^ ^d^d^ 


(6.16) 


Equations  (6.8)  and  (6.16)  yield 
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(6.17) 


or 


(6.18) 


The  noise  equivalent  resistance  of  the  device  is  defined  as 


(6.19) 


From  equation  (6.18)  one  obtains 


Equation  (6.20)  will  be  used  to  calculate  the  equivalent  resistance  as  a 
function  of  the  bias  voltage.  Before  presenting  the  experimental 
results  in  Chapter  VII,  the  equipment  used  will  be  briefly  discussed  in 
the  next  section. 


6. 2. a The  Spectrum  Analyzer 

The  spectrum  analyzer  used  was  the  HP  853A/8558A  system.  The  HP 
853A  is  the  spectrum  analyzer  display  and  the  HP  8559A  is  the  spectrum 
analyzer  itself.  The  HP8559A  spectrum  analyzer  employs  harmonic  mixing 
to  cover  a measurement  range  of  lOMHz  to  21GHz  in  six  frequency  bands. 


6.2  Description  of  the  Equipment 
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It  can  display  frequency  spans  as  narrow  as  lOOKHz,  and  as  wide  as 
9GHz.  The  HP  8559A  can  be  used  to  measure  signals  over  an  amplitude 
range  of  -llldBm  to  +30dBm. 

6.2.b  The  High-Frequency  Preamplifier 

The  high-frequency  preamplifier  used  in  this  experiment  was  a two- 
stage  cascode-pair  JFET.  This  configuration  provides  good  bandwidth 
with  reasonable  gain.  The  bandwidth  is  enhanced  by  the  use  of  3uH  in 
series  with  the  load  resistance.  It  has  a gain  of  30dB  and  a -3db 
frequency  of  30MHz.  This  amplifier  was  designed  by  C.F.  Whiteside 
[9].  The  circuit  diagram  is  shown  in  Figure  6.5.  A radio  frequency 
amplifier  model  QB-258  by  Q-Bit  Corporation  was  cascaded  to  the  previous 
ampl i fi er . 

In  Chapter  VII  the  experimental  results  will  be  presented. 


54 


CHAPTER  VII 
EXPERIMENTAL  RESULTS 


Measurements  were  performed  on  the*l.lum  n-type  device  at  room 
temperature  at  11  and  17MHz.  It  was  pulsed  biased  up  to  IV.  The  device 
resistance  measured  by  noise  measurement  as  described  earlier,  gave  a 
value  of  r^j  = (.65  ± .06n),  which  is  very  close  to  the  DC  measurement 
performed  during  the  current-voltage  characteristic  measurement. 

Equation  (6.20)  from  Chapter  VI  is  used  to  calculate  the  noise 
equivalent  resistance  R^,  of  the  device.  In  Figure  7.1  R^  is  plotted  as 
a function  of  the  applied  voltage  for  Tg  = 300K  and  fg  = 11.0  MHz  and  fg 
= 17.0  MHz.  For  bias  voltages  up  to  .2V  the  measured  value  of  R^  is 
equal  to  the  small  signal  resistance  of  the  device  r^,  within  the 
experimental  error.  For  bias  voltages  above  .2V,  R^  increases  sharply 
with  the  voltage  applied.  These  results  will  be  compared  to  numerical 
simulation  using  the  Monte  Carlo  method.  This  will  be  the  subject  of 
the  next  two  chapters. 
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Figure  7.1:  Noise  equivalent  resistance  as  a 

function  of  the  applied  voltage 


CHAPTER  VIII 
MONTE  CARLO  SIMULATION 

8.1  Introduction 

In  order  to  calculate  transport  coefficients  such  as  electron 
mobility  for  a given  applied  external  field  one  has  to  determine  the 
distribution  function.  The  distribution  function  changes  in  response  to 
external  driving  fields,  the  perturbed  distribution  function  is  governed 
by  the  appropriate  Boltzmann  equation,  and  hence  solution  of  this 
Boltzmann  equation  provides  a theoretical  description  of  the  behavior  of 
the  system.  Ideally,  the  Boltzmann  equation  can  be  solved  by  analytical 
means;  and  then  one  has  an  explicit  mathematical  formula  for  properties 
of  the  electron  system  in  terms  of  the  governing  parameters  and 
variables:  for  example,  mobility  versus  magnetic  field  as  a function  of 
effective  masses  and  electron-lattice  coupling  constants  [16].  This  is 
possible,  however,  only  for  some  particular  cases  of  high  symmetry  and 
simplicity,  which  seldom  occur  in  practice. 

Simulation  techniques  have  been  developed  in  theoretical  physics 
and  engineering  to  avoid  these  difficulties.  The  Monte  Carlo  method  is 
one  of  these  simulation  techniques  which  proved  to  be  immensely  powerful 
in  electron  transport  theory. 

Basically,  a Monte  Carlo  transport  simulation  consists  of  modeling 
the  motion  of  an  electron  as  a sequence  of  free  flights  (in  an  applied 
field)  interrupted  by  collisions  (or  scatterings).  Between  collisions, 
the  electron  obeys  classical  laws  of  motion  determined  from  the  band 


57 


58 


structure  (effective  mass  and  nonparabol i city)  of  the  material.  The 
collisions  are  regarded  as  random  events  whose  probabilities  are  known 
functions  of  energy.  The  duration  of  a free  flight,  the  kind  of 
scattering  process,  and  the  change  in  momentum  produced  by  a scattering 
are  determined  using  numbers  produced  by  a pseudo  random  number 
generator  usually  available  on  the  computer  itself.  In  effect  this 
method  produces  a numerical  solution  to  the  Boltzmann  transport  equation 
for  a specified  set  of  scattering  processes  and  material  parameters 
[17].  Figure  8.1  shows  the  main  steps  in  a Monte  Carlo  simulation:  A 
carrier  is  initialized  by  giving  a momentum  Pq  with  a lattice 
temperature  T.  The  carrier  then  accelerates  freely  under  an  applied 
field  F until  it  is  determined  that  an  interaction  has  occurred;  this 
determination  is  based  on  total  scattering  rate  as  a function  of  carrier 
energy  and  on  a pseudo  random  number  selected  at  the  beginning  of  the 
free  flight.  After  a free  flight,  another  pseudo  random  number  is 
selected;  this  number  and  the  individual  scattering  rates  determine 
which  of  various  interactions  occurs.  Once  a scattering  mechanism  has 
been  selected,  additional  pseudo  random  numbers  determine  the  state  of 
the  carrier  for  the  next  free  flight.  These  steps  (free  flight  ended  by 
a scattering  event)  are  repeated  until  a specified  number  N of 
interactions  has  occurred  or  until  a specified  time  has  elapsed.  In 
fact,  the  history  of  an  electron  is  followed  for  a very  long  time  T and 
each  parameter  is  averaged  in  time.  In  the  steady  state  case,  this  time 
average  is  equal  to  the  ensemble  average.  Each  of  these  steps  will  be 
discussed  more  in  detail  in  the  following  paragraphs. 
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Figure  8.1:  Flow  chart  for  Monte  Carlo  simulation 
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8.2  General  Method 

8. 2. a Generation  of  Random  Numbers  with  Given  Distributions 

Usually  parameters  like  free  flight  time  between  collisions  are 
random  variables  whose  probability  distribution  functions  are  known.  It 
is  convenient  to  express  those  random  variables  as  a function  of  a 
pseudo  random  number  provided  by  the  computer.  These  random  numbers 
have  usually  a uniform  distribution  between  0 and  1.  The  transformation 
goes  as  follows: 

Let  r be  the  pseudo  random  number  and  f(r)  its  probability  density 
function.  In  this  case  f(r)  = 1 for  r E [0,  1]  and  f(r)  = 0 elsewhere. 

Let  0 be  the  random  variable  representing  a certain  physical 
quantity  and  p(9)  its  probability  density  function.  We  map  the  pseudo 
random  number  r to  the  distribution  function  of  the  physical  quantity  0 
as  follows: 


P{r'  < r}  = P{0  < 0^}  (8.1) 

or 

Iq  f(r‘)dr'  = /q^  p(0)d0  (8.2) 


Since  f(r)  =1  r E [0,  1] 

9 

r = Iq^  p(9)d9 


(8.3) 


From  equation  (8.3),  9^  can  be  expressed  as  a function  of  the 
pseudo  random  number  r.  This  method  is  used  to  generate  quantities  such 
as  free  flight  time,  scattering  angle.  The  transformation  resulting 
from  equation  (8.3)  is  sketched  in  Figure  8.2. 
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8.2.b  Initial  Conditions 

In  the  steady  state  situation,  the  time  of  simulation  must  be  long 
enough  so  that  the  initial  conditions  of  the  electron  motion  do  not 
influence  the  final  results.  The  choice  of  a "good"  time  of  simulation 
is  a compromise  between  the  need  for  ergodicity  and  the  request  to  save 
computer  time  [18].  If  a highly  improbable  value  of  the  wave  vector  k 
of  the  electron  is  chosen,  the  first  part  of  the  simulation  can  be 
strongly  influenced  by  this  unappropriate  choice. 

A common  practice  is  to  choose  as  initial  energy  kgT  where  kg  is 
the  Boltzmann  constant  and  T the  lattice  temperature  and  deduce  from 
that  the  initial  wave  vector  kg  as 

I^qI  =1  (2m*kgT)^/2 


8.2.C  Electron  Motion 

In  a simple  semiconductor  electrons  or  holes  are  considered  as  free 
particles  with  effective  mass  m*  and  a dispersion  law 

2i*-  = Y(e)  (8.5) 

In  the  case  of  a parabolic  band  which  is  considered  here  y(e)  = e where 
e is  the  energy  of  the  carrier. 

The  effective  mass  m*  is  usually  considerably  smaller  than  the 
free-electron  mass  mg  = 9.1  10"^^  kg.  In  the  case  of  GaAs  it  is 
necessary  to  describe  the  electron  with  two  different  effective  masses. 
Between  collisions  the  electron  obeys  the  law  of  motion 
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Figure  8.2:  Generation  of  stochastic  variable 

9 with  a given  distribution  function 
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- ef  (8.6) 

where  -e  is  the  electron  charge,  f the  electric  field.  If  the 
electric  field  P is  parallel  with  the  k^-axis,  we  have  the  situation 
shown  in  Figure  8.3. 

If  At  is  the  drifting  time,  equation  (8.6)  can  be  rewritten  as 

k^(t)  = k^.  - At  (8.7) 

where  is  the  value  at  the  start  of  the  electron  flight  and  k^  is 

the  value  at  the  end  of  the  flight. 

Since  the  k-space  of  the  electron  is  symmetric  about  the  electric 
field  ? , it  is  useful  to  have  an  equivalent  two-dimensional 

(kp,  k^)  -space  where  kp  = (k^  + instead  of  the  three- 

dimensional  (kj(,  ky,  k^)  [19].  Hence  the  initial  state  would  be  (kp, 
k^J  and  the  final  state  (kp,  k^^)  and  only  the  z-component  of  equation 
(8.7)  needs  to  be  considered.  So 

kz^  = ; (8.8) 

since  ? = (0,  0,  -F)  , each  free  flight  is  terminated  by  a 

scattering.  Since  there  are  several  competing  mechanisms  the 
appropriate  scattering  mechanism  is  chosen  by  another  Monte  Carlo 
procedure.  The  free  flight  duration  At  has  to  be  chosen  by  a Monte 
Carlo  procedure  too.  These  will  be  discussed  in  the  next  paragraph. 
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Figure  8.3;  Electron  motion  in  momentum  space 
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3.2.d  night  Duration  and  Choice  of  Scattering  Mechanism 

Each  scattering  mechanism  is  characterized  by  a transition  rate 
i<' ) from  the  state  l<  to  the  state  . The  subscript  n 
denotes  an  individual  scattering  process  and  can  take  values  n = 1,  2, 
N if  there  are  N mechanisms  considered.  The  total  scattering  rate 
from  the  state  i<,  due  to  the  n-th  process,  is 

\i^)  = / di<'  (8.9) 


So,  the  total  scattering  rate,  due  to  all  processes,  is 


^ N 

X(k)  = E X (i<)  (8.10) 

n=l  ^ 


where  i<  is  a function  of  time.  Usually  the  total  scattering  rates  are 
only  functions  of  k = |i<|  , so  x(i<)  is  easily  transformed  to  x(E). 

Let  assume  that  the  electron  drifts  for  a time  t before  being 
scattered  and  that  this  time  consists  of  n tiny  increments 

Atp  At2,  At^.  The  probability  of  the  electron  being  scattered, 
within  the  time  interval  At-  , is  X(i<)  At-  , where  x(i<)  is  the  total 
scattering  rate  defined  by  equation  (8.10).  The  probability  that  there 
will  not  be  any  scattering  during  At^  is  (1  - x(j<)  At^).  Since 

those  are  independent  events,  the  overall  probability  of  an  electron  not 
having  any  scattering  during  a time  t is  then 


n 

P(t)  = n (1  - x(i<)  At.)  (8.11) 

i=l  ^ 


In  {P(t)}  = E In  {(1  - x(ic)  At.)} 
i=l  ^ 


(8.12) 
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Since  x(j<)  At.  <<  1 


In  {(1  - x(ic)  At^)}  ~ - x(i<)  At^ 


(8.13) 


and 


n 


In  {P(t)}  = - Z X(j<)  At. 

i=l  ^ 


(8.14) 


or  finally 


P(t)  = exp  (-  /q  X(j<)  dt') 


(8.15) 


The  probability  density  of  an  electron  drifting  for  a time  t and  then 
suffering  a collision  is 


Using  equation  (8.3)  which  maps  the  probability  density  p(t)  to  a 
uniformly  distributed  random  number  r we  obtain 


Mathematically  it  is  possible  to  invert  equation  (8.17)  and  express  t as 
a function  of  the  random  number  r.  In  general  it  is  complicated  to 
perform  such  an  inversion  since  an  integral  equation  needs  to  be  solved 
for  each  scattering  event.  The  difficulty  is  not  fundamental  but  lies 
in  the  numerical  evaluation  of  the  integrals.  Rees  [20]  has  devised  a 
very  simple  method  to  overcome  this  difficulty.  He  introduced  a new 
fictitious  "self-scattering"  such  that  the  total  scattering  probability, 
including  this  self-scattering,  is  constant  and  equal  to  a certain 
constant  P . 


p(t)  = X(i()  exp  (-  /q  X(i<)  dt' ) 


(8.16) 


r = 1 - exp  (-  /J  x(i<)  dt') 


(8.17) 
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r = x(i<)  = Xj^(ic)  + x^{t)  (8.18) 

where  Xj^(i<)  is  the  total  real  scattering  rate  and  X^(j<)  the  total 
virtual  or  self-scattering  rate.  If  the  carrier  undergoes  such  a self- 
scattering, its  state  after  the  collision  is  taken  to  be  equal  to 

its  state  k before  the  collision  so  that,  in  practice,  the  electron  path 
continues  unperturbed  as  if  no  scattering  at  all  had  occurred. 

Generally  it  is  enough  that  ^ is  larger  than  the  maximum  of  x(t)  . 

N 

Here  r is  taken  as  equal  to  Z where  the  X were  defined 

n inax  n 

in  equation  (8.9). 

If  A(i<)  is  replaced  by  P in  equation  (8.17),  then 


r = 1 - exp  (-P^) 

(8.19) 

1 - r = exp  (-rt) 

(8.20) 

t = -iln  (1  - r) 

(8.21) 

However,  since  r is  evenly  distributed  between  0 and  1,  so  also  is  (1  - 
r),  and  in  practice  the  following  expression  is  used: 

t = - In  (r)  (8.22) 

The  time  given  by  equation  8.22  is  not  usually  the  flight  time  between 
real  scattering  processes,  but  is  more  likely  to  be  the  time  between  one 
virtual  scattering  process  and  another.  The  computer  time  "wasted"  in 
taking  care  of  self-scattering  events  is  more  than  compensated  for  by 
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the  simplification  of  the  calculation  of  the  free-flight  time 

duration.  The  free  flight  must  be  terminated  by  a real  scattering 

event.  A further  random  number  is  then  used  to  select  a mechanism  among 

all  n different  possible  scattering  mechanisms.  This  random  number  r is 

i 

compared  to  the  partial  sums  E normalized  by  P , where  i 

n=l 

varies  from  1 to  N the  total  number  of  scattering  mechanisms  considered 

i .. 

in  the  model.  If  r is  smaller  than  1 E then  the  i^"  scattering 

mechanism  is  selected.  If  r is  greater  than  the  total  scattering  rate 
N 

iE  then  a self-scattering  or  virtual  scattering  is  chosen. 

n=l 

This  scattering  process  selection  is  described  in  Figure  8.4.  In  fact 
the  method  is  exactly  similar  to  that  described  earlier  by  equation 
(8.3).  Instead  of  a continuous  probability  distribution  function,  we 
have  here  discrete  events  which  are  described  in  terms  of  sum  of  6 
functions.  The  integral  should  be  replaced  by  a discrete  sum  and  the 
equality  by  a inequality  sign. 

8.2.e  Choice  of  the  State  after  Scattering 

Once  a type  of  scattering  is  selected  to  end  a free-flight  a new 
state  i<  must  be  chosen.  If  the  free  flight  ended  with  a virtual 
scattering,  the  state  must  remain  the  same  before  and  after  that  virtual 
scattering.  One  the  other  hand,  if  a real  scattering  occurred,  then  the 
state  i<  must  be  generated  stochastically  according  to  the  differential 
cross  section  of  that  particular  mechanism  [18].  Especially  different 
techniques  are  used  to  choose  the  post  scattering  angles  for  randomizing 
and  non-randomizing  mechanisms.  The  scattering  involving  acoustical 
phonons  is  randomizing  while  those  involving  polar  optical  phonons  are 
not.  These  techniques  will  be  described  in  more  detail  in  section  8.3 
when  the  Monte  Carlo  method  is  applied  to  calculate  transport 
coefficients  of  a particular  semiconductor  material:  GaAs. 
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N 


Self  -scattering 


Figure  8.4:  Scattering  channel  selection  with 

N real  scattering  processes 
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8.2.f  Collection  of  Results 

The  data  collected  at  each  free  flight  will  form  the  base  for  the 
determination  of  the  quantity  of  interest. 

We  obtain  the  average  value  of  a quantity  A[k(t)]  (e.g.,  the  drift 
velocity,  the  mean  energy,  and  so  on)  during  a single  history  of 
duration  T as 


<A>t  = Y /J  A[k(t)]dt  (8.23) 

or 

1 

<A>t  = T ^ ^ A[k(t)]dt  (8.24) 

where  the  total  flight  duration  time  T has  been  split  into  free  flights 
of  duration  t^ . When  a steady  state  quantity  is  computed,  T should  be 
taken  sufficiently  long  enough  so  that  <A>j  in  equation  (8.24) 
represents  the  ensemble  average  of  the  quantity  according  to  the  ergodic 
theorem. 


8.3  Application  to  Bulk  GaAs  Semiconductor 

The  energy-momentum  band  structure  (E(k)  versus  ± k)  of  the  GaAs 
material  is  shown  in  Figure  8.5.  There  is  one  central  valley,  with  a 
small  radius  of  curvature,  where  electrons  have  a low  effective  mass 
equal  to  0.067  mg,  and  there  is  a shallow  satellite  valley  at  the 
equivalent  (100)  Brillouin  Zone  Edges.  By  symmetry,  however,  identical 
sattelite  valleys  are  located  at  the  (010)  and  (001)  edges  as  well.  The 
valence  band  is  well  separated  from  the  conduction  band  so  the  transport 
of  holes  is  not  considered. 
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For  every  type  of  scattering  mechanism,  the  post  scattering  energy 
E(k')  is  calculated  such  that  the  energy  conservation  law  is  observed. 
In  case  of  scatterings  where  acoustic  phonons  are  involved,  the  process 
is  assumed  to  be  elastic  with  E(k)  = E(k'). 

The  list  of  scattering  processes  considered  in  the  work  is  given  in 
Table  8.1  with  their  respective  transition  rate  Sp(k,k'),  total  scatter- 
ing rate  ^r|(E)  3nd  the  energy  conservation  law  for  each  type 

of  scattering  mechanism.  The  material  parameters  are  given  in  Table  8.2. 

8. 3. a Angle  Selections 

As  mentioned  in  section  8.2.e  the  state  of  the  electron  after  a 
scattering  process  depends  on  the  type  of  process  involved.  For  the 
case  of  acoustic  and  intervalley  processes,  the  transition  probability 
Sp  (k.k‘)  (as  listed  in  Table  8.1)  does  not  depend  upon  the  angle 
between  k and  k'  so  that  all  wave  vectors  that  satisfy  the  appropriate 
law  of  conservation  of  energy  are  equally  probable. 

If  the  final  momentum  state  k'  makes  an  angle  0 with  the  k^-axis 
direction,  then  the  probability  density  p(9)  is  proportional  to  the 
number  of  states  on  the  circumference  of  a circle  of  radius 

Ik*  I sin  9 = (2m*  E'  sin  0 (8.25) 

where  E‘  is  the  energy  of  the  electron  after  scattering.  This  is 
sketched  in  Figure  8.6. 

We  assume  that 

p(9)  = A sin  0 (8.26) 


The  normalization  requires  that  A = 1/2. 
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Figure  8.5:  Band  structure  of  GaAs 
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Table  8.1:  Scattering  mechanisms  in  GaAs  from  Boardman  [19] 


Scattering  Transition  Rate  Total  Scattering  Rate 

Mechanism 


(h/2i)S(K,K') 


X(k)  e X(E) 


Acoustic  phonon  for  hD^  | N 6(A) 

either  absorption  or  5 

emission  (intervalley)  (2os)(2r)'' 


2 4 

4nD  s^h 


^A  ~ hs|K-K' I 


absorption  or  emission 
A = E'-E  = 0 


Polar  optical 
phonon  { i ntraval 1 ey ) 


2fip^hu) 
0 


X 

4nKjj(2n)^ 


/ 2h(4rK^)^  E^^^  jjl/2  - ^,1/2 


X = absorption 

= ( Hg *1 ) 4 ( Ae ) emission 

Ng  = [exp(hv,^/RgT)-l]-^ 


y 


= absorption 
= (^g'^1)  emission 
= E'  - E ± h<v^  = 0 


Equivalent  intervalley 
( satel 1 i te-satel 1 i te) 

(z-l)hD^  X 

[z-i).fhy/h 

(2x)^(2pu,g) 

3 

/ 2 upw  h 

X = ^g4(^g)  absorption 

y = Ng  absorption 

= i5(Ag)  emission  (N  +1) 
Ng  = [exp(hWg/k|^T)-l]" 

= emission  (N^+l) 

A = E*  - E t hu) 

a ,e  e 

Non-equivalent  interval  ley 

ahD^x 

3m*  3/2g2g,l/2Y 

a)  central  satellite 

20Wg(2m)^ 

X = N^5(A^)  absorption 
= (lig'''l)4(Ag)  emission 

Y 

''a  e 

/ 2tipu 

P 

= absorption 

= (N^+1)  emission 
* E'  - E + 4 ± hiv^ 

b)  satel 1 i te  - central 

hO^X 

n 

m*3/2g2g,l/2, 
c n 

2oiv^(2n)^ 

/ 2Trpu) 

X = absorption 

Y 

= N^  absorption 

(Ng+l)6(Ag)  emission 

(N^+l)  emission 

Ng  = [exp(h(v^/kgT)'h'^ 

^.e 

= E'  - £ - 4 t huj 

1 
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Table  8.2:  Data  for  GaAs  from  Boardman  [19] 


Density,  p 

Velocity  of  sound,  s 

High  frequency  dielectric  constant, 

Static  dielectric  constant, 

Polar  optical  phonon  frequency, 

Equivalent  intervalley  phonon  frequency, 
Non-equivalent  intervalley  phonon  frequency, 

Acoustic  deformation  potential  in  central 
and  satellite  valleys,  D^ 

Equivalent  intervalley  deformation  potential,  D^ 
Non-equivalent  intervalley  deformation  potential,  D^ 
Central  valley  effective  mass,  m* 

Satellite  valley  effective  mass,  m| 

Valley  separation,  a 


5.37  g/cm^ 

5.22  X 10^  cm/s 
10.82 
12.53 

5.37  X 10^^  rad.s"^ 
4.54  X 10^^  rad.s"^ 
4.54  X 10^^  rad.s"^ 

7eV 

lO^eV/cm 
lO^eV/cm 
0.067  mQ 
0.35  mQ 
0.36  pV 
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Mapped  with  a uniform  random  distribution,  sin  6 can  be  generated 
using  the  following  relations: 


M 

0 CD 
X) 
CD 

Q. 

CD 

(8.27) 

r = 1/2  (1  - cos  0) 

(8.28) 

the  new  state  k'  is  therefore 


i<'  = ( |k'  I sin  0,  |k'  I cos  9) 


(8.29) 


The  situation  for  polar  optical  phonon  scattering  is  completely 
different,  since  the  transition  rate  Sf^(k,k')  as  shown  in  Table  8.1 
depends  upon  the  angle  between  k and  k'  through  the  factor  |i<  - . 

It  is  necessary  to  consider  first  the  angle  3 between  k and  k' 
using  k as  reference  as  shown  Figure  8. 7. a. 

The  probability  density  p(S),  as  in  the  previous  mechanisms  is 
proportional  to  sin  3 and  also  must  be  weighted  by  the  factor  |k-k'|“^. 


p(3)  = 


A sin  3 
|k  - 


(8.30) 


O 

The  factor  |k  - k'|  can  be  expressed  in  terms  of  E and  E',  the 
energies  of  the  states  k and  k'  respectively  and  3 as  follows: 


|k  - k'|2  = (E  + E'  - 2(EE‘)^/2  3) 


then 
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l**x 


Figure  8.6:  Angle  selection  for  acoustic  and 

intervalley  process 
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Figure  8.7:  Angle  selection  for  polar  optical 

phonon  process 
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P(e)  = 


8 sin  S 


E + E'  - 2(EE')^''^  cos  3 


where  B = 


A 

2^ 


Let  f = 


2(EE' 


(e1/2  . e-1/2)2  ’ 


X = cos  3 


Then  cos  3 can  be  chosen  by  a random  number  r as 


r = 


rl  dx 

•*cos  3 (l  + f - fx) 
(-1  dx 

^-1  (1  + f - fx) 


or  carrying  out  the  integrals 


- (1  + f - f cos  3) 

In  ( 1 + 2f) 


Thus 


cos  3 = Ul/.  ^ 


(8.32) 


(8.33) 


(8.34) 


(8.35) 


As  is  shown  in  Figure  8. 7. a,  all  possible  positions  of  ]<'  are 
generators  of  a cone  about  i<  , of  semi-angle  3.  The  azimuthal  angle  <<) 
shown  in  Figure  8. 7. a can  be  selected  with  a random  number  such  that 


= 2tt  r (8.36) 

r E [0,  1] 

The  coordinates  of  it'  in  the  - kp  plane  is  then 

k'^  = |k' I cos  Y (8.37) 
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k' p = |k' I sin  Y 

where  y is  the  angle  between  and  f as  shown  in  Figure  8.7.b.  The 
angle  y is  related  to  ()>  and  8 in  the  following  way: 


cos  Y = cos  8 cos  + sin  8 sin  ij;  cos  (|>  (8.38) 

where  ip  is  the  angle  between  1<  and  ? . The  angle  t|)  can  be  determined 
by  cos  ^ = k^/|k|  where  is  known  at  the  end  of  the  free  flight  and 
|k|  = . 

For  the  different  types  of  scattering  considered  here  then,  the 
post-scattering  wave  vector  !<'  is  completely  determined.  The  wave 
vector  i<'  will  be  the  initial  value  for  the  next  free  flight. 

8.3.b  Drift  Velocity  and  Diffusion  Coefficient 

The  velocity  of  an  electron  can  be  expressed  as 


so  the  mean  velocity  is 


<v> 


fi/m* 


/an  i:  p ^ i< 
/all  f P (^)  <1  i< 


(8.40) 


However,  in  the  k^  direction  the  motion  is  uniform  since  the  electron  is 
under  a constant  electric  field  f , then  the  distribution  p(k^)  = 1 and 
if  the  integral  is  split  into  a sum  of  integrals  over  a free  flight 
equation  (8.40)  becomes 
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<v> 


1i/ni*  E ^ 

free  flights 


d k 


z 


(8.41) 


This  represents  the  mean  drift  velocity  in  the  z-di recti  on.  Carrying 
out  the  integrals,  we  obtain; 

''drift  ^ 2m*  0 (8.42) 

'^Zf 

where  D = E /.  d k 

^ . 

This  is  very  convenient  to  use  because  k_  , k^  are  z-components  of  the 

ZLf 

final  and  initial  k vector  for  a free  flight.  They  are  known  once  the 
energy  is  known  and  the  type  of  scattering  mechanism  is  selected.  The 
computation  of  v is  then  a matter  of  recording  these  values.  The  drift 
velocity  as  a function  of  electric  field  obtained  from  this  method  for  a 
bulk  GaAs  semiconductor  is  shown  in  Figure  8.8. 

The  longitudinal  diffusion  can  be  obtained  from  Einstein's  equation 


var  z = 2 D t (8.43) 

where  z is  the  displacement  and  D the  diffusion  coefficient. 

It  is  therefore  expected  that  if  the  motion  of  a single  electron 
from  an  initial  position  (0,0,0)  is  observed  over  a sampling  time  T and 
a large  number  of  samples  are  taken,  then  the  displacements  for  the 
different  samples  should  have  the  same  variance. 
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To  make  good  use  of  the  simulation  of  the  drift  velocity,  equation 
8.43  can  be  rewritten  for  a fixed  T as 

T var  = 2 D (8.44) 

or 

D = 0.5  T (<v^>  - <v^>^)  (8.45) 

where  is  the  average  drift  velocity  in  the  z direction  and  T sampling 
time. 

The  simulation  of  the  drift  velocity  can  be  performed  in  two 
ways:  either  the  number  of  real  scattering  events  is  fixed  or  the  total 

sampling  time  is  fixed.  The  latter  method  is  used  to  calculate  the 
diffusion  coefficient,  since  to  be  valid,  equation  (8.43)  requires  that 
T is  constant.  The  diffusion  coefficient  D as  a function  of  the 
electric  fields  is  plotted  in  Figure  8.9  for  a bulk  GaAs  semiconductor. 
The  programs  used  to  simulate  v(E)  and  D(E)  are  listed  in  Appendix 

1. 

The  drift  velocity  and  the  diffusion  coefficient  obtained  from  this 
Monte  Carlo  simulation  are  used  to  calculate  the  high  frequency  noise. 
This  will  be  described  in  Chapter  IX. 


o Monte  Carlo  simulation 
^ Exper i menta 1 data 


■J 1 1 1 I I 1 I L. 

1 23456789 


Figure  8.8:  Electron  drift  velocity  as  a function 

of  electric  field 
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Figure  8.9:  Electron  diffusion  coefficient  as  a 

function  of  the  electric  field 


CHAPTER  IX 
DISCUSSIONS 

In  Chapter  VII,  the  experimental  results  of  the  high  frequency 
noise  measurements  were  presented.  In  Chapter  VIII  the  electron  drift 
velocity  and  the  electron  diffusion  coefficient  were  computed  using  the 
Monte  Carlo  simulation  method.  One  would  like  to  compare  the 
experimental  data  with  the  theoretical  calculations.  That  will  be  the 
subject  of  this  chapter. 

Using  Poisson's  equation  and  the  current  equation,  neglecting  the 
displacement  current  and  the  diffusion  current,  one  can  express  the 
length  of  the  device  L and  the  applied  voltage  at  the  anode  Vg  as  a 
function  of  the  applied  field,  in  an  integral  form; 


„ _ eA  r^L  E^v(E)dE 

0 " I ^0  AN„Eu(E) 

° (1-q-^i ) 

0 


(9.1) 


L 


Eu(E)dE 
AN.Ey(E) 
(1-q-^T ) 


(9.2) 


where  Vg  = bias  voltage,  Ig  = DC  current,  A = cross-sectional  area,  L = 
contact  spacing,  E|_  = electric  field  strength  at  the  anode,  Nq  = 
concentration  of  fully  ionized  donors,  and  e = permittivity.  In  these 
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two  equations  the  field  dependence  of  the  electron  mobility  is  taken 
into  account.  Moreover,  Gisolf  and  Zijlstra  [21]  derived  an  integral 
form  of  the  noise  spectral  density,  taking  space  charge  injection  and 
carrier  heating  into  account; 


D(E)(EL-E)^dE 
ANny(E)E  . 
(1-q-T 


(9.3) 


where  D(E)  is  the  electron  diffusion  coefficient.  Equation  (9.3) 
together  with  equations  (9.1)  and  (9.2)  forms  a system  of  integral 
equations  with  three  unknowns,  namely  E[_,  u(E),  and  D(E).  This  system 
is  difficult  to  solve  analytically.  Numerical  techniques  are  required 
to  circumvent  such  difficulty.  This  will  be  discussed  in  the  following 
section. 


9.1.  Numerical  Decomposition  of  the  Current-Voltage 
Characteristic  and  the  Noise  Spectrum  Density 

Bosman  [22]  developed  a method  to  solve  numerically  equations 
(9.1),  (9.2),  and  (9.3).  He  considered  equally  spaced  bias  voltage 
values  V„.  (i  = 1...N).  For  each  V„  , there  is  a one  to  one 
correspondence  between  the  current  value  Iq.,  the  noise  spectral  density 
SAVq^  , and  the  electric  field  strength  E|  . . To  every  E[_.  corresponds  a 
value  of  u(E)  and  D(E). 

He  assumed  that  between  two  consecutive  values  of  Ei  , the  electron 
mobility  and  the  electron  diffusion  coefficient  are  constant; 


D(E)  = D^-  and  (E)  = 


(9.4) 
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for  E,  < E < E, 

i ^i+1 

Using  the  assumptions  stated  in  equation  (9.4),  equations  (9.1),  (9.2), 
and  (9.3)  can  be  rewritten  as 


eA 


i 

z 

s=l 


UsE^dE 

(l=q-4^) 

°i 


L = 


eA 


i 

E 

S = 1 


UgEdE 

AfOTJ- 

(1-q-^) 

°i 


(9.5) 


(9.6) 


S,Vi(f) 


i£9A  ^ 

? F 

^0.  ^=1  L 

1 


D^(E^_-E)‘'dE 

^ 1 M 3 

S-1  (1-q-j ) 


°i 


(9.7) 


The  integral  equations  are  now  transformed  into  a system  of  three  non- 
linear equations  with  three  unknowns,  E|_.,  and  . For  a given  i, 

equations  (9.5)  and  (9.6)  are  solved  first  to  obtain  and  Ei  , then 

■-i 

these  values  are  substituted  in  equation  (9.7)  to  calculate  . All  the 
integrals  in  equations  (9.5),  (9.6),  and  (9.7)  can  be  evaluated 
analytically.  This  process  is  done  iteratively  to  give  u^,  , and  E^. 

in  ascending  order  of  i.  The  details  of  these  calculations  are  given  in 
Appendix  2. 

The  accuracy  of  the  results  depends  strongly  on  the  magnitude  of 

the  interval  [E.  , Ei  ].  Another  difficulty  in  the  calculations  is 

^i+1 
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the  fact  that  the  numerical  solutions  can  be  unstable  due  to  the 

presence  of  a pole  in  the  integrand  of  the  integrals  in  equation  (9.5), 
(9.6),  and  (9.7). 

The  diffusion  coefficient  and  the  electron  drift  velocity  obtained 

from  this  decomposition  are  plotted  along  with  the  values  obtained  by 

Monte  Carlo  simulation  in  Chapter  VIII.  A difficulty  one  might 

encounter  in  the  measurement  of  high  field  properties  of  GaAs  is  the 

possible  occurrence  of  traveling  domains.  These  domains  would  make  the 

electric  field  inside  the  device  dependent  on  position  and  time. 

However,  according  to  Kroemer  [23],  these  dipole  domains  are  formed  only 

when  Np  L exceeds  2 x 10^2  cm“2.  This  is  not  the  case  here  where  Nq  is 
15  -3 

2 X 10  cm  and  L is  1.1pm.  No  traveling  dipole  effects  were  observed 
during  the  measurements. 

In  conclusion  of  these  high  frequency  noise  measurements  we  found 
that  the  noise  of  a 1.1pm  n-type  GaAs  device  at  high  bias  voltages  can 
be  explained  in  terms  of  the  hot  electron  properties  of  the  mobility  and 
diffusion  coefficient  if  the  field  dependence  of  these  parameters  is 
taken  into  account  properly. 

In  Chapter  X,  we  will  present  the  conclusion  of  the  overall  work 
and  suggest  some  numerical  simulations  and  experimental  work  to  pursue 
in  these  very  short  structures. 


CHAPTER  X 

CONCLUSION  AND  RECOMMENDATIONS  FOR  FUTURE  WORK 
10.1  1/f  Noise  and  Ballistic  Transport 

10.1.  a Numerical  Simulation 

One  problem  encountered  in  Chapter  V for  the  calculation  of  was 
the  determination  of  the  total  number  carriers  Ngf^r.  As  we  should 
understand,  ^eff  total  number  of  carriers  which  suffer 

scatterings.  For  short  devices,  as  Hesto  [24]  mentioned,  an  important 
fraction  of  the  electrons  does  not  suffer  any  collision  at  all.  For  a 
given  electric  field,  this  fraction  of  ballistic  electrons  can  be 
obtained  through  Monte  Carlo  simulation.  Hesto  showed  at  77  K,  this 

fraction  is  still  important  for  distances  as  long  as  0.35  ym.  This 
consideration  would  improve  greatly  the  computation  of  Ng^^  mentioned  in 
Chapter  V.  Such  a Monte  Carlo  simulation  may  give  information  about  the 
current  due  to  ballistic  and  non-bal 1 istic  electrons.  Device 

simulations  initiated  by  Tehrani  [14]  should  be  further  pursued  taking 
into  account  these  ballistic  effects.  The  idea  is  as  follows; 

The  Monte  Carlo  simulation  should  give  th  electron  distribution  at 
each  time,  then  the  electric  field  is  calculated  by  solving  Poisson's 
equation. 

10.1. b  Other  Experiments  on  Short  Structures 

As  stated  previously  in  section  10.1. a,  for  low  temperature  and  if 
the  thickness  of  the  low  doped  region  is  sufficiently  small 
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(L  < 0.5  m),  an  important  fraction  of  the  electron  population  can 
flow  through  the  device  having  a small  number  of  interactions  or  even  no 
interaction  at  all  with  the  polar  optical  phonons  and  then  can  reach  the 
anode  with  high  kinetic  energy.  It  has  also  been  shown  in  Chapter  II 
that  the  current-voltage  characteristic  of  a device  is  not  sufficient  to 
determine  whether  the  transport  is  ballistic  or  not. 

An  experiment  which  allows  separation  of  groups  of  electrons 
according  to  their  kinetic  energy  is  probably  interesting. 

Hesto  et  al . [25]  proposed  the  following  experiment.  They  added  a 
Schottky  contact  at  the  end  of  the  regular  N'*’  N“  N'*'  device  in  a way  that 
it  opposes  its  variable  potential  barrier  <))g  to  the  electron  energy. 
The  device  is  described  in  Figure  10.1. 

Let  respectively  Vg^  and  V^g  be  the  N“  and  the  Schottky  bias 
voltages.  The  electrons  having  a ballistic  flight  through  the  N"  region 
have  a fast  transit  through  the  N"*"  region  and  then  move  against  the 
retarding  electric  field  of  the  Schottky  barrier.  By  this  way,  only 
those  having  their  kinetic  energy  at  the  end  of  the  N region  higher  than 
the  Schottky  barrier  4>g  reach  the  anode.  By  varying  V^g  one  controls 
the  height  of  this  Schottky  barrier.  The  result  would  result  in  good 
separation  of  these  groups  of  electrons. 

10.2  High  Frequency  Noise  Measurements 

The  same  measurements  as  described  in  Chapter  VI  and  VII  should  be 
repeated  for  the  shorter  devices.  However,  the  simulation  of  v(E)  and 
0(E)  should  be  changed.  The  method  used  in  Chapter  VIII  for  Monte  Carlo 
simulation  is  not  valid  anymore  in  the  case  of  short  devices.  The 
ergodicity  theorem  does  not  apply  in  this  case.  Instead  of  time 
average,  we  have  to  perform  an  ensemble  average.  The  average  value  is 
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Figure  10.1: 
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Schottky  barrier  device 
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performed  at  a preselected  time  so  that  we  would  have  the  average 
velocity,  the  diffusion  coefficient  as  a function  of  time. 

Taking  a Fourier  transform  of  the  diffusion  equation  Jacoboni  and 
Reggiani  [18]  showed  that  the  diffusion  coefficient  can  be  calculated 
through  Monte  Carlo  simulation  as  a function  of  frequency  and  the 
reciprocal  of  the  distance.  This  approach  is  necessary  if  one  wants  to 
calculate  the  diffusion  coefficient  at  short  times  and  short  distances. 


APPENDIX  1 

MONTE  CARLO  SIMULATION  PROGRAM 


ooc~>  ooo 


COMPUTES  SIMULATION  OP  THE  0PI=T  V E L C I T ' Y , M = 4 N 
ENERGY  AND  DI3TRI5UTI0N  =UNCTI0N  Or  1-0  ' 

ELECT  ROM 

D0U3.L;  ?R£CiSICN*6  X?  ,HRHC/K2i=  /KZI/KT  ,M 
1/fiO/NE/Nl/MTOY,VSVG/V24VC/’V3AVG/V4AVG,CIr 
1/  Ir4,=T/TF/7I,TT0T 

DCU3LE  ?R:CISI0K*E  P ( 1 C ) , G A MM 4 ( 2 ) , E = ( 1 0 ) , SC  A TT ( 1 0 ) 

1  /OK  ( 2) /L (1 C) 

INTEGER'S  M = ^VV/VR/TT,V/SR/SS/</GY,  4Y,NCI  = 

PUNOAMENTAL  CONSTArjTS 

H=1.05459 
E = 1 . 6021  R 
C = 4. 3C321  3 
K3=1.3£C62 
M=9. 10913O 

OATA  ON  MATERIAL  GAAS 
WRITE(6/41  ) 

41  PORMATdOX/'MATERAL  0 EM  S IT  r ( GM  . C M- 3 ) ' ) 

READ (1 1 / 639) RhE 

WRITE  (6/42) 

42  PORMATCIOX/ 'VELOCITY  OP  SOL ND  ( 1 0* » 3 C M , S- 1 ) ' ) 

READ  (1 1 ,639) s 

WRITE ( 6/ 43) 

43  format (lOX, 'HIGH  FREOUENCY  DIELECTRIC  CONSTANT') 

READ ( 1 1 ,639) R1 

WRITE(6,44) 

44  POkMAT(10x,'L0W  FRcDuemcy  CIELECTRIC  CONSTANT') 

REA J ( 1 1 , 639  ) R 2 

WRITE(6,1  ) 

1 P OR.M  AT  ( 1 OX  , 'OP‘^1  C AL  =HONQf|  P R E 0 U E N C Y ( 1 0 * * 1 3 R A 0 . S ~ 1 ) ' ) 
KEiG(11/o£9)W0 

WfsiT 

2 P0R.MAT(10X,'EDU.  INTERVALLEY  PHONON  PREG.') 

R£AD(11,639);jc 

WRITE (5,3) 

3 FOR  MATdOX, 'IN  TER  VALLEY  PHCNON  FRECUENCY') 
READ(11,639)wi 

W R I T E ( 5 , 4 ) 

4 FORMAT ( 1 CX, ' ACOUS TIC  DEFORMATICN  ? C T E NT  A L ( = V ) ' ) 

READd  1 , 633)  THA 

WRITE (6, 5) 

5 FORM  AT  ( nx,  ' EGU  . INT  ERV  alley  0 E F . oqT  E N T I A L ( E V ) ' ) 

READ ( 1 1 ,639) THE 

W R I T E ( 6 , 6 ) 

6 format(igx,'intekvalley  oef. potential') 

REA0(11,6S9)THI 
W R I T E ( 6 , 7 ) 

7 FOR  MATdOX, 'VALLEY  SEFAPATIOM') 

REAOd  1,639)0 

WRITE(6,3) 

3 F0R-1AT(10X, 'CENTRAL  VALLEY  EFFECTIVE  MASS') 
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READ (1 1 , 639) 3M1 
WRIT£(o/9) 

9 FOR, MATC10X, 'SATELLITE  VALLEY  EFF.'iasS') 

REAO ( 1 1 ,639) EMI 

639  FORMATCEIO. A) 

C 

C 

C FINAL  OATi  input 

r 

Ao  WRITE(5,1  D) 

10  FCRiMAT(10X,'TEMP£RAT'JR£(KELVIN)') 

REAOC11,531)  T 

5 31  F 0 R M A T ( £ 1 0 . A ) 

WRITE (5,21 ) 

21  FOR MATClOX, 'MAXIMUM  ENERGY(EV)') 

RcAO(11,532)  cMax 

532  FORMAT(EIO.A) 

WRI TE ( 5,  1 1 ) 

11  FORMATdOX, 'NUM3ER  OFAVERAGES') 

18  R£AO(11,533)  NCIF 

533  FCRMATCI5) 

IFCTHI.NE.O.O)  go  to  750 
IF(iOIF.LE.  300C)  go  to  16 
WRITE(5,  710 

710  FCRMATCIOX,' MAXIMUM  COLLISIONS  IN  CNE  vaLL=Y=TC0O') 

GO  TO  1 S 

750  IF C NOIF. lE. 20020)  GO  TO  16 
W K I T E < 5 , 1 7 ) 

17  FORM  A r C 1 Ox, ' M A X IMUM  COLLISIONS  IN  TWO  V A L L E Y S = ’ C 00  ' ) 

GO  TO  13 

16  WRITE ( 5, 1 2 ) 

12  F0RMAT(10X,'nUM3£R  OF  ELECTRIC  FIELOS') 

R £ A 0 ( 1 1 , 5 3 4 ) N F 

5 3 4 F 0 R • A T ( I 2 ) 

00  30  1=1, NF 
WRITE(5,13)  I 

13  F0R.MAT(10X,'  = :EL0',I2) 

READ  ( 1 1 , 5 3 5 ) ( I) 

5 35  FOR.MATC  El  2.  A) 

30  CONTINUE 

W R I T E ( 5 , 2 3 ) 

23  FORMA T ( 1 Ox, ' V ALL E Y FOR  0ISTRI3UTI0N  FUNCTION  :1  OR  2') 
RcAD(11,536)  VV 

536  F0RMAT(I2) 

WRITE (5, 22) 

22  form  AT ( 1 0X, '01 STA NCE  FR0«  <2  AXIS  IF  0ISTRI3UTI0N  :1  TO  21' 
REA0(11,537)  VR 

537  FCRMAT(I2) 

REA3(11,532)  OL 

C 

C CALCULATE  PHONON  FPEOUENCIES  AND  OCCUPATION  RATIOS 

H',iC=H*WO/(E-*100.0) 

hwi=h*wi/(e*ic:.0) 

HwE=H*wE/ (E*10C.0) 
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IFIWO.NE.0.3)  GO  TO  909 

N0=0.0 

GO  TO  903 

909  NO=l/(EX?(C^O*76.335)/T)-1) 

908  NI=1/(£X?((WI*76.385)/T)-1) 

NE=1/(£XP((W£*76.335)/T)-1) 

C 

C CONSTANTS  P09  SCATTERING  RATES 

C 

C1=1.0£-^12*C*C*SQRT(M)*WC*(1/R1-1/K2)*(N0  + 1)/(1.4U2*H*SQRT(£)) 
C2=C1*N0/ (NO+1 ) 

C3=1.3E+10*(2*,M)**1.5*K3*T*THA*THA*E*E*SCRT(e)/(4.0*3.142*RH5 
1 •S*S'»H*H*H*H) 

C4=2.0£+14*M**1.5*THE*THE*E*E*(ME+1)*SCRT(E) 

1/(1.4U2*3.142*RHE*WE*H*M*h) 

C5=C4*NE/ (NE+1 ) 

Co=1.0£+14*(EMl*M)**1.5*THI*THI*E*e*(NI+1)*SCRT(E)/(1. 4142*3. 142 

1 «RHc*'iJl*H*H*H) 

C7=3.0£+14*(EN2*M)**1.5*thI*THI*E*E*(NI+1)*SQRT(E)/(1.4142 

1*3.142*RH£*WI*H*H*H) 

C6=C7*NI/ (NI+1 ) 

C9=C6*NI / ( NI+1 ) 

C 

C CALCULATE  X SPACE  i*ESH  ELEMENT  FOR  90H  VALLEYS 

C 

0K(1)  = 1.0£  + 7*S':RT(2*EMl*M*eMAX*E)/(H*20.0) 

Dn(2)  = 1.0E  + 7*S':pT(2*EM2*M*EMAX*E)/(H*2C.O) 

C 

C CALCULATE  VALUES  OF  ,<  Z AT  CENTRE  OF  EACH  MESH  ELEMENT 

C IN  CHOSEN  valley  AND  WRITE  TO  FILE 

C 

XZI=-0K(VV)*9.3 
00  305  LL=1/20 
WRI TE ( 5/ 807)  K ZI 
807  FOR.MATCEI  2.4) 

KZI=KZI+DX(VV) 

305  CONTINUE 

C SET  PARAMETERS  >=0R  CENTRAL  VALLEY/THEN  CALCULATE  THE 

C TOTAL  SCATTERING  RATS  FCR  REAL  FROCESSES(R)  FOR  A NUM3ER 

C OF  ENERGIES  UP  TO  THE  MESH  SIZE. STORE  MAXIMUM  VALUE  OF  R 

C IN  GAMMA(I)  TO  CALCULATE  PSEUDO  (SELF)  SCATTERING  RATE. 

C TT=1  ENA3LES  PROGRAM  TO  RETURN  TO  LA3EL  40 

C 

TT=1 
EM=EM1 
V = 1 

31  GAMMA (V) =0. 0 

£1=0.0 

J = 1 

35  EI=EI+EMAX/20. 0 
GO  TO  100 

40  R=0.0 

CO  50  1=1/10 

R=R*L(I) 

50  CONTINUE 
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c 

c 

c 


71 

75 


72 

C 

C 

C 

c 


30 


u 

C 

c 

c 


c 

c 

c 


IrCR.OT.GiMMaCV) ) GAMMa(V)=R 
J = J + 1 

IFCJ.NS.21)  GO  TO  35 

SET  PiRA meters  for  SATELLITE  VALLEY  AND  REPEAT  PROCESS  TO 
05TAI.M  GAMMA  (2)  . 

IrCV.Ea.E)  GO  TO  71 
EM  = EM  2 
V = 2 

GO  TO  31 

WRITEC5/75)  GAMMA(1)/GAMMA<2) 

WRITE(7/75)  GAMMa(1)/GAMMA(2) 
FORMA7(10Xx'GAMMa(1)=',c10.4/'GAMMA(2)='/=10.4) 

WRITE  < 5/ 72 ) T, NO  IF, EMAX 
WRITEC7/72)  T,NOIF,£MAX 

FORMAT(1GX/'TE.MP='/£10.4/5X,'REALCLSNS  = '/  I4/5X/'MAX=NERGY  = ' 
1/E10.4) 

SET  MESH  REGISTERS  TO  ZERO  AND  PLACE  ELECTRON  AT 
STARTING  POINT  IN  M£SH.TT=0  FOR  ITERATIVE  PROCESS 


TT  = 0 
J = 1 

TTOT=O.C 
0T0T=0.0 
NN  = 0 

VAVG=0.0 

V2AVG=0.0 

0IF=0.0 

V3AVG=0. 0 

V4AVG=0.0 

0IF4=0.Q 

CriECX=0.0 

TF=4.0E-10 

FT=TF 

I0IF=0 

ET0T=0.0 

XZTOT=O.Q 

N P R I N T = 0 

V = 1 

EM=EM1 
XRH0  = 0.  j 
SR  = G 
SS  = 0 
GMAX=0 

IF  NO.  OF  REAL  PROCESSES  EQUALS  CHOSEN  VALUES 
ENO  ITERATION  ANO  GO  TO  FINAL  CALCULATION. 


CALL  RANDOM  NUMB E R < NO T =C ) ANO  CALCULATE  TIME  OF 
FLIGHT  UNOEP  ELECTRIC  FIELD  ANO  NEW  POSITION  OF 
ELECTRON  X SPACE 
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37 
9 0 


C 

c 


94 
C 

C 

C 

C 

C 

95 


96 


C 

c 

c 

c 

r 

C 

W 

c 

c 


1 00 


R=9ANU(0.0/1 .0) 

KZF=(1Q.O/H)»SOR7(2Q.O*K3*T*5M*M*(-4LOG(R))) 
go  to  a/o 

R=RANU (0.0/1 .0) 

IF(R. £0.0.0)  ?=1.0E-00 
TIM£=-AL0G ( R ) / GAMMA ( V ) 
k:i=kzf-k(tim  = *=,  = CJ)  *1 .0£  + 1 3)  / H 
KT=S3RT(KRH0*KRH0+KZI*<ZI) 

£I=r(*H»KT*KT*1.0£-14/(=*2*£H*:M) 

IF  electron  leaves  mesh  place  it  on  edge  of  mesh  ANO 
REGISTcR  OCCURRENCE  IN  COUNTER  GMAX 

IF(£I.LE.£MAX)  GO  TO  95 
GMAX=GMAX+1 
E I = £ M A X 

KT=1.QE-*-7*S3RTC2*EM*H*EMAX*E)/H 
IF(KZI.GT.O.O)  GO  TO  94 
KZI=-SaRT(A3S(XT*KT-KRH0*XRH0)) 

GO  TO  95 

KZI=30RT(AES(KT*XT-KRHO*KRH0)) 

STORE  FLIGHT  TIME  IN  TOTAL  TIME  REGISTER  FQR  APPROPRIATE 
VALLEY/THEN  REGISTER  PASSAGE  OF  ELECTRON  THROUGH  ELEMENTS 
OF  X SPACE  MESH 

ttct  = ttot+-time 

IF(NN.NE.O)  GO  TO  96 
0 = H » 1 0 . 0 » T ri  E • ( K Z I + X I F ) / ( 2 * E M * M ) 

0TOT=0T0T+0 

IFCOTOT. LT.OL)  go  TO  96 
tf=ttot-time 

FT  = TF 
NN=NN+1 

XZT0T  = XZT0T  + ASS(X2I-XZ'=) 

£TOT=ETOTr(x:i*x;i-XZP»XZF)»10.0*H/(2»EM*M) 

SUM  TOTAL  changes  IN  X SPACE  POSITION  ANO  ENERGY  SPACE  POSITON 
ANO  STORE  IN  X2T0T  ANO  ETCT.SUM  MEAN  ENERGY  CHANGE  IN  MNE. 

SUM  INOIVIOUAL  -values  FOR  EACH  VALLEY  IN  XZ(V)  AND  ET(V). 


CHECX  FOR  ROUNDING  ERRORS  LEADING  TO  NEGATIVE  ENERGY  VALUES. 
1=  this  OCCURS/PLACE  ELECTRON  AT  STARTING  POSITION. 


IF(EI.GT.O.O)  GO  TO  100 

KRHO=0. 0 

KZF  = 1 . Gc  + 6 

£M=EM1 

V = 1 

GO  TO  90 

CALCULATE  FINAL  ENERGY  VALUE  FOR  EACH  SCATTERING  PROCESS 


£F(  1 ) =EI-HNO 
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EF(2)=£I+HW0 
EF(3) =£1 
EF(4) =EI 
£F(5.)  = EI-rtWE 
£F(6)  = £I'*Tt^(£ 

EF(7)=£1-MWI+0 
EF(a)=£I-HWI-0 
EF(9) =EI+HWI-D 
EFC10)=EI+HWI+0 
C 

C ' SCATTERING  RATES  FOR  REAL  PROCESSES 

C 

IF(EF(1  ) .GT.-G.O)  go  to  110 

L(i)=o.a 

GO  TO  120 

110  L(1)=C1*SQRT(EM)*AL0G(ABS((SCRT(EI)+SQRT(EF(1)))/(SQRT(En 
1SQRT(EF(1)))))/SQRT(£I) 

C 

C EMISSION  OF  OPTICAL  PHONON 

C 

120  IF(£F(2) .GT.0.0)  GO  TO  125 

L(2)=0.0 
GO  TO  130 

125  L(2)=C2*SORT(EM)*ALOG<ABS((SQRT(£I)+SQRT(5P(2)))/(SORT(EI) 
1 SORT (EF ( 2) n ) ) /SORT  (El) 

c 

C ABSORPTION  OF  OPTICAL  PHONON 

C 

130  IF(EF(3) .GT.O. C)  GO  TO  135 
L(3)=0.0 
GO  TO  140 

135  L(3)  =C3*£M**1  . 5*S3RT(E'=(3)  ) 

C 

C EMISSION  OF  ACOUSTIC  PHONON 

C 

140  IF(£F(4) .GT.G.C)  GO  TO  145 

L(4) =0.0 
GO  TO  150 
145  L(4)=L(3) 

C 

C ABSORPTION  OF  ACOUSTIC  PHONON 

C 

150  IF(EF(5) .GT.0.0)  GO  TO  155 
L(5)=0.0 
GO  TO  170 

155  IF(V.£Q.2)  GO  TO  160 

L(5)=0.0 
GO  TO  170 

160  L(5)=C4*EM**i,5*sqrt(EF(5)) 

C 

C EMISSION  OF  EOUIVALENT  INTERVALLEY  PHONON 

C 

170  IF(EF(6) .GT.0.0)  GO  TO  175 
L (6) =0. 0 
GO  TO  150 
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175 

IFCV.50.2)  GO  TC  >:0 
L(6) =0.0 
GO  TO  190 

180 

c 

L(6).=C5*“M**1  .5*S:9T(zF(6)) 

C 

c 

ASSOnFTIO'J  of  SvUIVaLENT  INT59VALL5Y  “HONON 

1 90 

IF(Gr(7).GT.O.O)  GO  TO  195 
L (7) =0.0 
GO  TO  210 

195 

IrCV.cu.2)  GO  TO  200 

L(7)=0.0 

GO  TO  210 

200 

r 

L(7)=C6»SQ9T(EF(7)) 

C 

c 

EMISSION  OF  INTER  VALLEY  FHONONCSJTELLITE  TO 

CENTRAL) 

210 

IF(EF(8) .GT.0.0)  GO  TO  215 
L( 3) =0.0 
GO  TO  230 

215 

IF(V.EO.I)  GO  TO  220 

L(3)=0.0 

GO  TO  230 

220 

c 

L(3) =C7*S0RT(EF(3) ) 

c 

c 

EMISSION  OF  INTERVALLEY  PHONON  (CENTRAL  TO 

SATELLITE) 

230 

IF(EF(9) .GT.0.0)  GO  TO  235 
L <9) =0 .0 
GO  TO  250 

235 

IF(V.EO.I)  GO  TO  240 

L(9)=0.0 

GO  TO  250 

240 

c 

L(91 =C8*SQRT(£F(P) ) 

c 

c 

ABSORPTION  Or  INTERVALLEY  ? H CNO N ( C EN T R A L TO 

SATELLITE) 

250 

IF(ErdO)  .GT.0.0)  GO  TO  255 

L(10)=Q.0 

GO  TO  270 

255 

IF(V.EQ.2)  GO  TO  260 

L(10)=0.0 

GO  TO  270 

2 6 0 
r 

L(10)=C9*SGRT(EF(10)) 

V 

c 

c 

ABSORPTION  OF  INTERVALLEY  PHONON (S ATELL HE  TO  CENTRAL) 

270 

r 

IFCTT.EO.I)  GO  TC  40 

w 

c 

c 

CALCULATE  SUM  OF  REAL  PROCESS  SCATTERING  RATE 

S 

280 

SCATT(1)=L(1)/GAMMA(V) 

00  290  !C  = 2/10 

SCATT(Y)=SCATT(l<-1)fL(K)/GAMMA<V) 

290 

CONTINUE 

100 


c 

C CALL  RANDOM  N UM B E R . S E L E C T SCATTERING  CHANNEL 

C 


R = RA,NU(0.0/1  .0) 

Ir(R.LT.SCATT(D) 

GO 

TO 

30C 

IF(R.LT.SCATT(2)) 

GO 

TO 

31C 

IF(R.lT.SCATT(3) ) 

GO 

TO 

3 2C 

IFCR.LT. SCATTC4) ) 

GO 

TO 

33C 

IF(R.LT.SCATT(5)) 

GO 

TO 

34C 

IFCR.lt. SCATTCS)  ) 

GO 

TO 

35C 

IFCR.lt, SCATTC7)  ) 

GO 

TO 

36C 

IFCR.LT.SCATTC3)  ) 

GO 

TO 

37C 

IrCR.LT. 3CATTC9)  ) 

GO 

TO 

33C 

IFCR.LT.SCATTCIC)  ) 

30 

TO 

390 

GO  TO  400 
C 

C SET  ENERGY  AFTER  SCATTERING  PROCESS 

C 

300  ErIN=EF(1) 

GO  TO  420 
310  ErIN=EF(2) 

GO  TO  420 
320  ErIN=EF(3) 

GO  TO  410 
330  £FIN=5F(4) 

GO  TO  410 
340  £FIN=EF(5) 

GO  TO  410 
350  £FIN=EF(6) 

GO  TO  410 
360  EFIN=EF( 7) 

GO  TO  430 
370  EFIN=EF(3) 

GO  TO  430 
380  EFIN=EF(7) 

GO  TO  430 
390  EFIN=EF(10) 

GO  TO  430 
400  £FIN=EI 

GO  TO  450 
410  SR=SR+1 
C 

C REGISTER  REAL  COLLISION. CALCULATE  K SPACE  POSITION  AFTER 

C ACOUST  IC / I NT EP V ALLE Y OR  EQUIVALENT  INTERVALLEY  PHONON  SCATTERING 

R=RANU(0.0/1 .0) 

KT=1.0£+7»SQRT(2*EM*M*EFIN*E)/H 

KZF=KT*(1-2*R) 

KRHO  = KT*SQRT  (4*R*  (1-Rn 
GO  TO  460 
420  SR=SR+1 
C 

C REGISTER  REAL  COLLISION. CALCULATE  K SPACE  POSITION  AFTER 

C optical  PHONON  SCATTERING 

C 
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c 

c 

c 

i*3Q 


440 


450 

C 

C 

C 


C 

C 

c 

c 

c 

4 <)0 


c 

c 

c 

c 

c 

470 


471 


R=RANU(0.0,1 .0) 

U=«ANU(0.0/1 .0) 

PHI  = 2-*3.142*R 

EX=i*SQRT(=FlN*EI)/((SQRT(EI)-SQRT<5FIM))**2) 

5£TA=(((1+EX)-(1+2*£X)**U)/=X) 

RHO=(3£TA*K2I/KT-SORT(A3S(1-5£TA»3£TA))*KRHO/<T*COS(PHT)) 

KT=1.0£*7»SQRT<2*£M*M*=FIN*£)/H 

X2F=XT*RHO 

KRH0=XT*SaRT(A3S(1-RH0*RH0)) 

GO  TO  A 6 0 

CHANGE  valley  PARAMETERS  FOR  INTEfiVALLEY  PROCESSES. 

IF(V.£2.1)  GO  TO  440 

V = 1 

£M=£M1 
GO  TO  410 

V = 2 

£M=EM2 
GO  TO  410 
SS=SS+1 


REGISTER  SELF  SCATTERING  PROCESS. K SPACE  POSITION  UNCHANGEO 
XZF=XZI 

KT=S2RT(KZF*KZF+KRH0*KRH0) 

EFIN  = H*H*KT  *KT»1  .0E~1 4/ ( E*2*EM*M) 

CHECK  IF  ELECTRON  IS  SCATTERED  OUT  OF  MESH. IF  SO/REGISTER 
PROCESS  ON  COUNTER  GMAX,AND  PLACE  ELECTRON  ON  EDGE  OF  MESH 
LA3EL  90  REPEATS  ITERATIVE  PROCESS  STARTING  WITH  FREE 
ELECTRON  FLIGHT  UNDER  ELECTRIC  FIELD. 

IF(EFIN.LE. EMAX)  GO  TO  90 
GMAX  = G.MAXtI 

KT=1  .0£  + 7*SCRT (2*EM*M*£MAX*E) /H 
IF (K RmO . GT . KT ) KRHO=KT 
KZF=S2RT(A3S(KT*KT-KRH0*KRH0)) 

GO  TO  90 

FINAL  CALCULATION  OF  DRIFT  V E L 0 C I T Y- V E L , T I M E SPENT  IN  EACH 
VALLEY-VTI,M(V),MEAN  ENERGY -MNE 


IF(IOIF.NE.O)  go  to  471 

VAVG=0.0 

V2AVG=0. 0 

V3AVG=0.0 

V4AVG=0 . 0 

IDIF=IDIF+1 

NPRINT=NPRINT+1 

£TOT=  = TOT-('<;2I*KZI-KZP*KZF)*1Q.O*H/(2*EM*M) 
KZT0T=KZT0T-A3S(KZI-KZF) 

KZI  = KZFt(=*f(J)*1.0E  + 13/H)'*((FT-TTCT)+TIME) 
£T0T=ET0T+(KZI*KZI-KZF*KZF)*10.0*H/(2*EM*M) 
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K2TOT=KZTOT+A3S(K2I-KZF) 

VI=cTOT/KZTOT 

VAVG=(VAVG*(IDIF-1)+VI)/IDIF 
V2AVG=(V2AVG*(ICIF-1)+VI*VI)/IDIF 
V3AVG=(V3AVG*(ICIF-1)+VI*VI*VI)/IDIF 
VAAVG=(V4AVG*(I0IF-n+VI*VI*VI*VI)/IDIF 
0IF=0.5  *FT*(V2AVG-VAVG*VAVG) 

DIF4=SGRT(A3S(V4AVG-4*V3AVG*VAVG+6«V2AVG*VAVG**2-3*VAVG**4)/12 
1 . 0 ) * ft 

CHECK=V3AVG-3*V2AVG*VAVG+2*VAVG**3 

£TOT=Q.O 

KZT0T=0. 0 

TTOT=0.0 

WRITE(5/430)'=(J)/VAVG/'0IF/IDIF 
WRITE(7/430)F(J),VAVGy0IF/IDIF 
430  F0RMAT(3(E9.3/1X)/1X,I5  ) 

NPRINT=C 

IFCIDIF.lt. NOIF)  GO  TO  87 
I0IF=0 
J = J + 1 

IFC J.NE.MF+1 ) GO  TO  30 

STOP 

END 


APPENDIX  2 


NUMERICAL  DECOMPOSITION  OF  I-V  CHARACTERISTIC  AND 
THE  NOISE  SPECTRAL  DENSITY 


For  a given  i,  the  sums  in  equations  (9.5)  and  (9.6)  can  be  written 


as 
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Equations  (A.l)  and  (A. 2)  can  be  rewritten  as 


(A. 2) 
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where 
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(A. 6) 


Note  that  the  subscripts  for  El  have  been  shifted  to  avoid  a subscript  0 
for  the  computer. 

Combining  equations  (A. 3)  and  (A.4)  after  carrying  out  the 

integrals,  we  can  express  Ei  as  a function  of  Ui 

'-i+1  ^ 
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(A. 7) 


Now  we  substitute  E.  (ui) 

^i+1  ' 


into  equation  (A.4) 
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We  have  to  find  the  zero  of  the  following  function  to  find  y^- 


Ey.dE 


(A. 9) 


Carrying  out  the  integral  we  have 


F(y.)  = C2  - L + eA  [ 


(A. 10) 


A simple  bisection  method  would  have  given  the  zeros  of  this  function. 
However,  in  some  cases  the  function  F(y^)  does  not  change  sign  and  the 
bissection  inethod  becomes  useless.  This  happened  when  F(uj)  has  a 
multiple  root  and  the  curve  is  just  tangent  to  the  y^-axis.  In  this 
case  we  have  to  check  if  the  derivative  of  F(y^-)  changes  sign  and  if  the 
absolute  value  of  F(y^-)  is  close  enough  to  zero.  If  these  two 
conditions  are  met  then  we  consider  a possibility  of  having  a zero  in 
the  neighborhood  of  that  point. 

The  derivative  of  F(y^)  is  calculated  by  using  Leibniz  formula 


E=E 
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Once  the  and  are  found,  D^-  is  easily  calculated  from 

equation  (9.7) 


S 


where 
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